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Abstract We survey the current situation regarding chemical modeUing of the syn- 
thesis of molecules in the interstellar medium. The present state of knowledge concern- 
ing the rate coefficients and their uncertainties for the major gas-phase processes - 
ion-neutral reactions, neutral-neutral reactions, radiative association, and dissociative 
recombination - is reviewed. Emphasis is placed on those key reactions that have been 
identified, by sensitivity analyses, as 'crucial' in determining the predicted abundances 
of the species observed in the interstellar medium. These sensitivity analyses have been 
carried out for gas-phase models of three representative, molecule-rich, astronomical 
sources: the cold dense molecular clouds TMC-f and L134N, and the expanding cir- 
cumstellar envelope IRC -1-10216. Our review has led to the proposal of new values and 
uncertainties for the rate coefficients of many of the key reactions. The impact of these 
new data on the predicted abundances in TMC-1 and L134N is reported. Interstellar 
dust particles also influence the observed abundances of molecules in the interstellar 
medium. Their role is included in gas-grain, as distinct from gas-phase only, models. 
We review the methods for incorporating both accretion onto, and reactions on, the 
surfaces of grains in such models, as well as describing some recent experimental ef- 
forts to simulate and examine relevant processes in the laboratory. These efforts include 
experiments on the surface-catalysed recombination of hydrogen atoms, on chemical 
processing on and in the ices that are known to exist on the surface of interstellar 
grains, and on desorption processes, which may enable species formed on grains to 
return to the gas-phase. 

Keywords Astrochemistry • Reaction Rate Coefficients ■ Gas-Phase Chemistry • 
Grain-Surface Chemistry • Chemical Modelling • Uncertainty Propagation • Sensitivity 
Analysis 



1 Introduction 

One of the major challenges in the field of astrochemistry is to understand fully the 
synthetic mechanisms that are responsible for the gaseous molecules that have been 
observed in the interstellar medium (ISM). Table 1 lists the molecules that have been 
definitely identified in the ISM, principally by observing spectra of two kinds: (i) tran- 
sitions between rotational levels observed in emission using ground-based telescopes 
that are sensitive at millimetre wavelengths, and (ii) transitions between vibrational 
levels observed in absorption at infrared wavelengths. Because of absorption in the 
Earth's atmosphere, many of the latter observations have been made with instruments 
mounted on satellites, such as the Infrared Space Observatory (ISO) and the Spitzer 
Space Telescope, which are also able to detect ices such as water ice on cold interstellar 
dust particles. Recent discoveries from the Herschel satellite have also been included. A 
large fraction, but by no means all, of the molecules listed in Table[l]are detected in the 
cold cores of dense interstellar clouds (ISCs) where, despite their name, the molecular 
density is quite low by terrestrial standards. The dominant gas-phase molecule in these 
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sources is molecular hydrogen, with the second-most abundant molecule, CO, down in 
abundance by a factor of 10"*. The high relative concentration of molecular hydrogen is 
due both to an efficient formation mechanism from hydrogen atoms on interstellar dust 
particles and the fact that hydrogen is the most abundant element in the universe. 

As well as being used to identify molecules that are present in the ISM, these 
spectroscopic observations can be used to estimate the abundances of the observed 
molecules, usually expressed as abundances relative to that of molecular hydrogen. 
The analysis of the rotational emission spectra is non-trivial due to complications 
such as heterogeneity of physical conditions in the source, the possibility of optical 
thickness in the spectral lines, and the fact that the analysis leads directly only to a 
column density along the line of sight in the upper (emitting) rotational level ( |Herbst| 
and van Dishoeck||2009 1. Converting this column density to a fractional abundance is 



itself a non-trivial matter since it is difficult to study molecular hydrogen, because it is 
non-polar. The resulting fractional abundances are then typically uncertain to at least 
a factor of ~ 3. 

Just as in atmospheric chemistry, attempts to further improve our understanding of 
how interstellar molecules are created in the observed abundances (that are estimated) 
involve two other interconnected kinds of scientific endeavor: (i) the construction of 
chemical models, which utilise large networks of elementary chemical and some physical 
processes, and (ii) experimental and theoretical efforts to obtain the information about 
these processes which is required as input to the models. Of primary importance are 
the coefficients which describe the intrinsic rate of all the elementary processes in the 
model, as well as the products of these processes. To define a rate coefficient, consider 
the process in which C atoms and H2 molecules radiatively combine to yield CH2. The 
loss of carbon atoms can be expressed by the differential equation: 

-d[C]/dt = fcRA[Cl[H2] (1) 

Here, k^A is the rate coefficient for radiative association, and [C] and [H2] are the 
instantaneous concentrations of C atoms and H2 molecules. Concentrations are usu- 
ally given in cm""^. Consequently, the units of fcj^A a-nd of the rate coefficients for 
all "second-order" processes are cm^ s~^. The values of rate coefficients generally 
depend on temperature and in astrochemical compilations it is usual to express the 
temperature-dependence of each rate coefficient by stating the values of three param- 
eters a, j3 and 7, where 

k{T) = a(r/300)'^ exp(-7/r) (2) 

However, only rarely is there sufficient experimental or theoretical information available 
to express the temperature-dependence of k{T) with much confidence. For instance, 
models often rely on extrapolations of k{T) out of its measurement range, which incurs 
very large uncertainties( Hibrard et al.|2009 l. 



The level of success of the models is tested by comparing the predicted molecular 
abundances with those derived from the spectroscopic observations. Despite some suc- 
cesses, the discrepancies that are found are frequently very large: greater than an order 
of magnitude. Of course, it must be recognized that both the predicted and observed 
abundances are uncertain. One recent development has been the attempt to quantify 
the size of the uncertainties in the predicted abundances by Monte Carlo methods, 
as discussed later in this article. Clearly, this procedure, along with estimates of the 
uncertainty in abundances derived from the spectroscopic observations, allows one to 
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decide whether or not a particular discrepancy between observed and predicted abun- 
dances is significant. A second development has been the use of models to identify 
'key' processes in the network; that is, elementary reactions which strongly influence 



the predicted abundances of one or more observed species and their precision ( Wakelam 



et al.|2005[|2006||Vasyunin et al.|2008||Wakelam et al.|2009[ ). At least, in principle, this 



information should encourage experimentalists and theoreticians to concentrate their 
efforts on reactions that are astrochemically crucial. 

The present review has its origins in a program supported by the International 
Space Science Institute in Bern, Switzerland. The authors of the current paper, along 
with some others whose contributions are acknowledged at the end of this article, 
held two meetings to discuss the current state-of-knowledge in respect of the syn- 
thesis of interstellar molecules. Some emphasis was placed on the importance of not 
only identifying the key astrochemical reactions, making the best possible estimates 
of their rate coefficients, but also of trying to quantify the uncertainties associated 
with these estimates so that some confidence could be placed in the significance of any 
disagreements between observed and predicted abundances. These efforts took place 
in parallel with a second related one: the development of a new database for astro- 
chemistry, named KIDA (Kinetic Database for Astrochemistry) which is now opened 
to the community : http://kida.obs.u-bordeauxl.fr/ Although estimates of uncertain- 
ties have been given in previous databases (e.g., Umist Database for Astrochemistry, 
UDfA http://www.udfa.net), KIDA goes further in that, for each key reaction, the 
entry includes a data sheet that gives reasons for the recommended value of the rate 
coefficient and its uncertainty, as well as references to any experimental and theoretical 
work on the reaction. 



2 History of chemical modeling 



Chemical modeling of sources in the ISM began in the 1950's following the detection 
of CH, CN and CH'*' via narrow absorption lines towards bright stars. In a landmark 



paper. Bates and Spitzer (19511 discussed quantitatively the density of molecules in 



interstellar space. Radiative association reactions, dissociative recombination (see Sec- 
tion [3|, and photodissociation and photoionization processes were emphasized, result- 
ing in a reasonable (if ultimately incorrect) explanation for the presence of CH+, but 
with questions concerning the formation of CH. About 10 years later, [McCrea and Mc-| 
Nally ( 1960 1 investigated the possibility of molecular formation by surface reactions 



on interstellar dust grains, focusing on H2, which had not yet been detected at that 
time, and CH. 

The detection of the A doubling transition of OH at 18 cm ( |Weinreb|1963| [Robinson] 
et al.|1964 l in the Galactic Center, followed by the detection at 6 cm of NH3, the first 
polyatomic molecule ( Cheung et al.|1968 l, and observational evidence for the presence 
of H2 towards e Per and ^ Per from rocket observations ( Carruthers|1970 1 stimulated 
interest in the molecular content of the interstellar medium and in the chemical and 
physical processes involved in the formation of the observed molecules. The advent of 
radio-millimeter telescopes and the subsequent discovery of many polyatomic species, 
including molecular ions, in the 1970's helped to put interstellar chemistry on a firmer 
basis. Not only astronomers but also chemists and physicists realized how promising 
and rewarding it could be to search for new, often exotic, molecules (see Table [T]) and 
to build chemical schemes with different formation and destruction processes. 
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Table 1 Reported Interstellar and Circumstcllar Gas-Phase Molecules on July 2010. 
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Initially, two main classes of regions were considered for quantitative chemical mod- 
eling: diffuse regions, with a total H density of ~ 10^ cm^'^ and a temperature of 50-100 
K, and dense regions, with an H2 density of ~ 10^ cm~'^ and a temperature of 10- 
20 K. The modeling was of a steady-state nature, in which the time derivative of all 
concentrations is set to zero. The detailed chemistry of diffuse regions, including the 
prototypical ( Oph cloud, was studied by Black and Dalgarno (19771, who included a 
detailed discussion of molecular hydrogen formation and excitation, and quantitative 
treatment of photoprocesses driving the chemistry. The chemistry of dense clouds, now 
known as cold cores, was investigated by Herbst and Klemperer (1973|, who showed 
that the chemistry is driven by cosmic rays, which initiate networks of ion-molecule re- 
actions. A major success was t he contribution of ch emical modeling to the identification 
of the IfCO^ molecular ion (Herbst and Klemperer 1974j . More recent observations 
have shown that both diffuse and dense sources are often parts of larger heterogeneous 
objects. Moreover, some cold dense cores are in the process of collapsing to form stars 
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and planetary systems. The collapse is a complicated one hydrodynamically, and can 
involve shock waves. 

The state of the field in the late 1970's was admirably described by two review pa- 
pers [DaIgam£an^^Iack]([T976| and I Watson] ( |1976| ); these describe the various physical 
and chemical processes relevant to reactions at the surface of interstellar grains and 
in the gas-phase and are still relevant today. A few years later, [Prasad and Huntress| 
(1980a I reported the first large-scale model for gas-phase chemistry, which was based 
on experimental studies of ion-molecule reactions performed with a variety of experi- 
mental techniques including their ion cyclotron resonance (ICR) facility. In addition, 
they provided an extensive list of rate coefficients for chemical reactions involving 
compounds of carbon, nitrogen and oxygen. Non-equilibrium effects were specifically 
addressed in a subsequent paper ( Prasad and Huntress|1980b 1. 

As far as interstellar modeling is concerned, one should also point to the decisive 
influence of computational power, which now makes it possible to solve a large num- 
ber of coupled chemical equations as a function of time from given initial conditions. 



In addition, numerical techniques, such as the Gear method (Gear 1971 Byrne and 



|Hindmarsh| 1987] ) , enable stiff problems, with time-dependent processes with orders of 
magnitude different time scales, to be handled. These techniques have been important 
for pseudo-time dependent chemical studies, which started with Leung et al. (19841. 
In this approach, the chemical evolution is solved independently of the physical condi- 
tions, which are taken as fixed and homogeneous, typically with an H2 density of 10* 
cm"'^ and a temperature of 10 K, conditions that are relevant to the "standard" TMCl 
cloud, often taken as a reference for the purpose of comparison between models and 
observations. Such numerical techniques are also used in the chemical interstellar shock 
models developed by Flower et al. (19851, who followed both the physical and chemi- 
cal evolution in a magneto-hydrodynamic (MHD) shock. Complex albeit steady-state 
models have been constructed for so-called Photon-dominated regions (PDR's), which 
are heterogeneous sources located near stars ( Hollenbach and Tielens|1997 RoUig et al 



20071. 



The number of species involved in such chemical models has increased to 400-500 
linked by over 4000 reactions, following the spectacular increase of detected interstel- 
lar and circumstellar molecules via their rotational (sub) millimeter- wave spectra. The 
need for chemical databases became urgent in the 1980's, and the group at Manchester 
( [Millar et al.|[T99T| ) was the first to publish such a database, which is now available 
on a website (UDfA at http: / / www.udfa.net /) , with regular updates published in the 
literature. Other databases are now also available including the OSU (Ohio State Uni- 
versity) database (e.g. http://www.physics.ohio-state.edu/~eric/research.html) and, 
as discussed later in this article, the already mentioned KIDA database. 

Consideration of both steady state and time-dependent studies has demonstrated 
the non-linear character of the chemical equations describing interstellar chemistry 
(Le Bourlot et al. 19931. Comparing steady state solutions with the integration of 



time-dependent equations has enabled the discovery of multiple solutions and bi-stable 
behavior. Two stable solutions have been shown to coexist within a typical range of 
control parameters, which can be " astrophysical" quantities (density, temperature, cos- 
mic ionization rate) or even the values of the reaction rate coefficients with a range of 
uncertainties. The solutions correspond to different chemistries adequately described 
by the level of ionization. The High Ionization Phase favors atomic ions, unsaturated 
species and radicals, whereas in the more standard Low Ionization Phase, small satu- 
rated molecules and molecular ions play a preponderant role. Such a possibility cannot 
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be excluded in the real ISM and it oflFers a natural way to explain current problems 
such as the low abundance of molecular oxygen as well as the amount of atomic carbon 
readily detected in dense clouds. 

A new stage was reached during the 1990's with the wealth of information obtained 
using the Infrared Space Observatory (ISO), launched by the European Space Agency 
in 1995 and operated for 3 years. Among the outstanding results obtained by the 
various instruments on board this satellite are the determination of the composition 
of the icy mantles of dust grains, the measurement of rotational transitions of H2 and 
other transitions such as the and O fine-structure lines, and the ubiquity of certain 
vibrational features thought to be due to polycyclic aromatic hydrocarbons (PAH's). 
These and other results have played a determining role in defining new needs for further 
developments in astrochemical models. 

The need for including surface processes, eventually triggered by photon irradiation, 
particle bombardment or thermal processing, has been raised, both to understand the 
composition of icy mantles ( Gibb et al.||2000 l and to find an alternative to gas-phase 
processes, which can fail to reconcile model predictions with observations for some 
specific molecules. An additional impetus for including gas-grain chemistry was the 
need to understand the chemical complexity resulting from the detection of increasingly 
complex organic molecules ( Hasegawa et aHl992) ). However, gas-grain models rely on 
many hypotheses concerning the energy barriers associated with both surface diffusion 
and surface reactions. These barriers are poorly constrained, and the energies that are 
assumed for them have a considerable impact on the results. An additional difficulty 



in gas-grain models arises from the stochastic nature of gas-grain chemistry (Tielens 



and Hagen||1982' I: under certain conditions, the rate of reaction is dominated by the 



accretion rate of species onto individual grains rather than by the rate of diffusion so 
that one should explicitly introduce the probability of the number of adsorbed species 
on the surface of one grain. Different ways have been considered to take these effects into 
account, including Monte Carlo techniques ( Charnley et al.|1997 Charnley|2001 1 and 
master equations ( Green et al.|2001 Biham et al.|20dT I . Lipshtat and Biham ( 2004b I 
and Barzel and Biham ( 2007 1 have shown subsequently that moment equations with an 
appropriate cut-off up to second-order could adequately describe gas-grain reactions. A 
simpler, but semi-empirical method has been introduced, in which the rate equations 
are modified to reproduce some of the effects of stochastic approaches ( [Caselli et al.| 
19981 [Garrodl|2008b| ). 



3 Gas-phase reaction rate coefficients and uncertainties 

Chemical models, designed to improve our understanding of how molecules are syn- 
thesised in different regions of the interstellar medium, consist of a large number of 
'rate equations'; that is, ordinary differential equations (ODEs), each of which, like 
eq. ([T|, represents how given chemical reactions or physico-chemical processes build- 
up or lessen the concentration of a particular chemical species. Currently, the OSU 
database osu-09-2008 includes 455 different species connected by 4457 rate equations. 
(The figures for the database for astrochemistry held at UDfA are very similar.) The 
455 species include 187 neutral atoms and molecules, 259 atomic and molecular cations, 
and 6 atomic and molecular anions. In addition, the concentrations of electrons and of 
neutral and negatively charged grains are included. In time-dependent models, assump- 
tions must be made concerning the starting conditions: the temperature, the elemental 
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Table 2 Number of reactions of different types in tfie gas-phase tliat are included in the OSU 
kinetic database for astrochemistry (osu-09-2008). 
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abundances and the state (atomic, molecular, ionic) of the elements, and then the 
ODEs are integrated forward in time to generate molecular abundances at different 
ages of the cloud. There are at least two requirements for any such model to be suc- 
cessful: (i) it should include all those processes that are important in determining the 
molecular abundances, and (ii) the values of the rate coefficients included in the model 
should be close to the correct values. It is the second of these two requirements that is 
considered, in general terms, in this section. 

Table [2] lists the different types of processes or reactions that are included in the 
OSU database and gives the number of each type of process in the model. The first 
type of reaction refers to the neutralization process between gas phase cations and 
negatively charged grains plus the H2 formation on grains. This latter process is ar- 
tificially included in gas phase models assuming that the H2 formation occurs at half 
the rate of hydrogen atom sticking on a grain ( Leung et al.|1984 |, i.e. the efficiency of 
the conversion is one as assumed by Jura (19741. Although the numbers of individual 
types of processes vary widely, it should not be inferred that these numbers reflect the 
importance of different types of processes. In this section, an attempt is made to sum- 
marize how the rate coefficients for some of the most important classes of reaction are 
estimated from the available experimental and theoretical evidence and, equally impor- 
tantly, how the uncertainties in the rate coefficients are estimated. We shall focus on 
the data for the following classes of reaction: (i) cation (and anion)-neutral reactions, 
(ii) neutral-neutral reactions, (iii) radiative association reactions (producing both ionic 
and neutral products), and (iv) dissociative recombination reactions. These categories 
include most of the reactions in astrochemical models. At the outset, it should be 
emphasized that the fraction of these processes for which rate coefficients 
have been measured or calculated at the low temperatures prevalent in cold 
cores is extremely small. In addition, for those reactions that may proceed to dif- 
ferent sets of products, the branching ratios to these different channels are frequently 
unmeasured. Consequently, the values of rate coefficients included in the models are ar- 
rived at by a combination of experiment, theory and common sense, although in many 
important cases neither experimental nor theoretical work has been performed, so that 
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chemical intuition, perhaps guided by information about chemically related reactions, 
may be the only way to estimate a rate coefficient. 

If the rate coefficient for a reaction is not to be negligibly small at the tempera- 
tures prevalent in cold cores, the value of 7 in eq. ([2| must be negative, zero, or at least 
very close to zero. In theoretical terms, this means that there should be no significant 
barrier on the path of minimum potential energy leading from reactants to products. 
In general terms, the value of the rate coefficient is likely to reflect the ability of the 
long-range attraction between the reactants to bring them into close contact. A sec- 
ond fundamental requirement is that there is an electronically allowed (or at least not 
strongly forbidden) path that leads from the reactants to products. Usually, this can 
be checked by applying the spin and orbital correlation rules to discover the symmetry 
of the potential energy surfaces that correlate with the separated reactants and the 
separated products ( |Smith|1980[ ). The importance of these rules is greatest when both 
reactants have one or more unpaired electrons. Then potential energy surfaces of dif- 
ferent spin multiplicity arise. The surfaces with the lowest spin multiplicity are likely 
to provide the lowest energy pathway from reactants to products. Collisions on other 
surfaces are unlikely to lead to reaction and this effect can introduce a factor < 1 into 
an estimate of the rate coefficient. 



3.1 Ion-Neutral Reactions 



The experimental database of rate coefficients for reactions between cations and neutral 



molecules is large. Following the early studies using ion cyclotron resonance ( Prasad and 



Huntress! 1980a Snow and Bierbaum|2008 1 and flowing afterglow methods, most of the 
recent kinetic data on ion-neutral reactions have been obtained using selected ion flow 
tubes (SIFTs) or some variant of this method. However, many of these measurements 
are only at room temperature or over a restricted range of temperatures around room 
temperature. The CRESU and ion-trap techniques have been applied to a limited 
number of reactions, yielding rate coefficients down to temperatures similar to those in 
cold cores. An important advantage of both these methods, as applied to ion-neutral 
reactions, is that mass spectrometry is used to observe changes in the concentrations 
of the reactant ions and the same technique can be used to identify the products of 
the reaction under investigation. 

However, for exothermic ion-molecule system, reactions between ions and polar 
neutral species have not been studied extensively in the laboratory. A few studies 
down to temperatures of ~ 20-30 K (Adams et al. 1985 Marquette et al. 1985 



Re- 



brion et al.|1988||Rowe|1988| ) and some studies at temperatures above 300 K ( |Viggiano| 
and Williams||200T I have been undertaken. Nevertheless, the rate coefficients of most 
thermal reactions used in model calculations must be estimated. These estimates can 
be used in the complete absence of experimental data or, for those cases where room 
temperature measurements are available, they can be used at higher and lower tem- 
peratures if there is reasonable agreement at 300 K. Detailed calculations of reaction 
rate coefficients require a knowledge of both the long-range and short-range aspects of 
the intermolecular potentials governing the reactions ( Clary|1988 1. For exothermic ion- 
molecule systems especially, it is often possible to neglect short-range considerations 
and to focus on the long-range asymptotic potential. 

The orbital angular momentum associated with the relative motion of the reactants 
under the influence of the long-range potential creates so-called centrifugal barriers. 
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Rather than solve the detailed quantum mechanics in the absence of short-range bar- 
riers, one can use a capture model, in which the sole criterion for reaction is that the 
reactants must pass over the centrifugal barrier to enter the short-range portion of the 
potential (Herbst 19961). Some additional dynamical barriers are created by anisotropics 



of the potential ( Troe||1987 l. The capture model thus provides an upper limit to the 



actual rate coefficient. For non-polar neutrals, the long-range potential depends solely 
on the distance R between ion and neutral, which in many instances is simply the R~'* 
potential arising from the interaction of charge and dipole polarizability. With this 
potential, the rate coefficient is given by the so-called Langevin expression: 



k-L = 27re(apoia/M)"^^^, 



(3) 



where e is the electronic charge, tjpoia is the polarizability, /i is the reduced mass of 
reactants, and cgs-esu units are utilized so that the units for the rate coefficient are 



For polar neutrals, the potential is anisotropic and depends on the orientation 
of the permanent dipole as well as the distance between reactants so that it is no 
longer central in nature. A number of detailed treatments of the long-range problem 
with the capture approximation have been undertaken, starting with the full quan- 
tum mechanical treatment. At a somewhat simpler level are the rotationally adiabatic 



approach of Clary (19871, the perturbed rotational states approximation of Sakimoto 



and Takayanagi ( |1980[ ), and the statistical adiabatic channel approach of Troe and 
co-workers ( Troe [l987| Maergoiz et al.|[2009 l. These detailed approaches can be used 
to calculate rate coefficients as functions of the rotational angular momentum of the 
neutral reactant or, by suitable averaging, thermal rate coefficients. 

A simpler approach for the astrochemist is based on classical trajectory calculations 
in which the long-range potential as a function of distance and orientation is used 
( Su and Chesnavich|1982 Maergoiz et al.|2009 l. Thermal rate coefficients determined 
from these calculations can be fit to analytical expressions depending solely on the 
reduced mass of the reactants, the polarizability, and the permanent dipole moment 
of the neutral reactant. This treatment is classical in nature and is only reliable at 
temperatures above which rotational motion can be considered classical, typically 10 
- 20 K. At lower temperatures, a semi-classical regime is entered, and at still lower 
temperatures 1 K), a fully quantum mechanical regime occurs, as studied in detail 
by Troe and co-workers (see, e.g., [Maergoiz et al.|2009| . The transition between classical 
and semi-classical behavior may be of some importance for cold cloud cores, where the 
temperatures lie in between the two regimes. The transition can be easily characterized 
(Troe 1989, 1996, Maergoiz et al. 2009) but, the classical expression is not seriously 
wrong until under 10 K. 

The best-known formula for the ion-dipolar rate coefficient /cq in the classical 



regime and for linear neutrals is the Su-Chesnavich expression (Su and Chesnavich 



1982 1. Maergoiz et al. ( 2009 and references therein) show that it works reasonably well 
for neutral species that are linear, symmetric tops, or asymmetric tops, although they 
have also explored the subtle differences among these structures. The Su-Chesnavich 



formula (see, e.g., Maergoiz et al.|2009 1 is based on the parameter x, defined by 



(SapolafeBT") 



1/2 ' 



(4) 



where /^d is an effective dipole moment of the neutral reactant, which is generally very 
close to the true dipole moment, and fc^ is the Boltzmann constant. The larger the 
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value of X, the larger the rate coefhcient. The expressions for /cq can be written in 
terms of the Langevin rate coefficient via 

ku/kL = 0.4767X + 0.6200, (5) 

if a; > 2, and 

/kh = (a; + 0.5090)^/10.526 + 0.9754, (6) 

if a; < 2. The latter reduces to the Langevin expression for x = 0. Alternatively, the 
expressions can be written in powers of temperature T. For example, for x > 2 

fco =ci+C2T-i/2 (7) 

where 



and 



(8) 



C2 = 2.1179-^i^. (9) 




For a; < 2, a more complex expression of the type 

ku=bi+b2T-^^^ + b3T-^ (10) 

pertains, where bi is simply the Langevin rate coefficient, while the parameters 62 and 
b3 are given by the expressions: 

62 = 0.4296-^^^ (11) 



and 

2 

^ 10.526fcB(«polaM)l/2' 

For a:: > 2, eq. (FtI) shows that if the second term is much greater than the first term, 

•—J —1/2 

which is more likely at low temperatures, the expression reduces to C2T ' , which is 
used in the OSU database for non-linear neutral reactants. For linear polar neutrals, 
the OSU database assumes that the species are rotationally relaxed (J ~ 0) and so uses 
the so-called "locked dipole" expression ( Herbst and Leung|1986 1 which is significantly 



larger than fe^. The UDfA database, on the other hand, uses the same temperature 
dependence but scales the result to a measured or estimated result at 300 K: 

fcD,UDfA = fcD(300 K)(300/T(K))i/2. (13) 

One should note that k^, does not diverge for T— )■ K as would be suggested by the clas- 
sical expressions (5), (7) and (13). Instead, levels off at a value fejj /ki^ = {l+fiY) /3apolaE 
where B is the rotational constant (in energy units) of the neutral ( |Troe|1987[ |1996| ). 



It is recommended that networks incorporate the Su-Chesnavich formulae discussed 
here in the situation where no experimental information is available, especially if infor- 
mation is required over a broad temperature range, say, from 10 - 1000 K. It must be 
recognized, however, that differences from the simpler T ' dependence may be small 
except for the prior use of the locked dipole approximation. This approximation should 
be removed from the OSU network for consistency, even if highly polar neutrals have 
their rotational levels relaxed thermally compared with the kinetic temperature. If a 
room temperature measurement is available, the capture theory result can be suitably 
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scaled. Use of the capture theory has been heretofore hindered by the lack of polariz- 
ability and dipole moment data, but a new compendium of information for interstellar 
neutral molecules calculated by Woon and Herbst ( 2009 1 makes the needed information 
readily available for all species in the networks, no matter how exotic. Instructions for 
determining the rate coefficients of ion-polar reactions in the range 10-1000 K can be 
found on the OSU URL. 

Most of the ion-neutral reactions that have been studied in the laboratory are 
between singly charged cations and neutral molecules that can be synthesised or pur- 
chased quite readily. It is clear from a study of Table [l] that these molecules comprise 
only a fraction of the neutral species that are potential reactants in dense clouds. 
The 'key' ion-neutral reactions discussed in Section |5.2| all involve radical atoms as 
the neutral reactant. To study these reactions involves an extra level of experimen- 
tal difficulty since the atoms need to be generated by flow-discharge techniques and 
their steady-state concentration estimated. The room temperature rate coefficients for 



a large number of such reactions have been recently tabulated (Snow and Bierbaum 



20081. 



A factor that might influence the rate coefficients of reactions between ions and open 
shell species (radicals) is the presence of spin-orbit (or fine structure) splitting in the 
case of radicals, which have one or more unpaired electrons. For example, the ground 
states of the C, N and O atoms are respectively '^P, *S and '^P. If the ion also has an open 
electronic shell, then potential energy surfaces (PES's) of different spin multiplicity will 
correlate with the reactants. This means that, in general, only a fraction of the collisions 
between the reactants will occur on a surface leading adiabatically to products. In these 
cases, it is recommended that a factor be introduced to allow for spin degeneracy. To 
take a specific example, the interaction between N(^S) atoms and C2H^ ions, which 
have a spin multiplicity of 2, will lead to PES's with both S = 1 and 5* = 2. It 
is reasonable to suppose that reactions will only occur on the triplet surface so any 
estimate of the rate coefficient should include a factor (3/8) to allow for this. Many 
radicals, atomic as well as molecular, possess electronic orbital angular momentum 
(i.e., L > 1 or yl > 1) and consequently the electronic state will be split by spin- 
orbit interaction into several components. However, as yet the influence of spin-orbit 
coupling on rate coefficients is not clear, we recommend that such effects be ignored. 



3.2 Neutral-Neutral Reactions 

It should be emphasized that, of the roughly 150 chemical species that have been 
identified in the interstellar medium (see Table [T]), the great majority are neutral and 
most of these are either unsaturated (hydrogen-poor) molecules or free radicals. The 
potential importance of neutral-neutral reactions in the chemistry of the interstellar 
medium, and in particular in cold cores, was underestimated for a number of years on 
the basis that many such reactions are known to have significant activation energies 
(that is, large positive values for 7 in eq. [2| with the result that their rate coefficients 
would be extremely small at the low temperatures encountered in the cold cores of 
dense clouds. Indeed, reactions between pairs of neutral molecules, e.g. H2 -f D2 — 
2 HD and C2H4 + C2H4 — >■ C4H8, are characterised by very high activation energies 
(corresponding to several thousands Kelvin) and can be ignored in models of most 
regions of the ISM. Many reactions between radicals and saturated molecules, e.g. D -|- 
H2 — >■ HD + H and CN -I- H2 — ^ HCN -I- H, have activation energies corresponding to 
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500 - 5000 K, and they can also be ignored in modeling the cold cores although there 
are a few exceptions; e.g., CN + NH3 HCN + NH2. 

All of the 382 neutral-neutral reactions included in the OSU database (osu-09- 
2008) involve, at least one unstable species; that is, one free radical. In only 99 of these 
reactions, is the other reactant a stable species: for example a saturated molecule, 
such as H2, CH4, or H2S, or an unsaturated molecule, such as C2H2 or C2H4. The 
reactants in the other 283 reactions are both unstable. The experimental database for 
the latter class of bimolecular reactions is much sparser than for the former, particu- 
larly at low temperatures, reflecting the experimental difficulties in measuring the rate 
coefficients for radical-radical reactions. For many reactions of both these types, a start- 
ing point for estimating the rate coefficients at low temperature can be the extensive 
reviews and evaluations of kinetic data for gas-phase reactions that have been com- 



piled for the purposes of modelling the terrestrial atmosphere (Atkinson et al. (2004 



2006 2007 2008 1; http://www.iupac-kinetic.ch.cam.ac.uk and Sander et al. (20061 



http:/ /jpldataeval.jpl. nasa.gov/download. html). 

Essentially all the kinetic studies of neutral-neutral reactions at temperatures below 
100 K employ either the continuous flow CRESU technique or the pulsed Laval nozzle 
variant of the technique. In a recent review ( Canosa et al.|20'08 1 , it is stated that 73 re- 
actions have been studied using the continuous CRESU method to reach temperatures 
as low as 13 K, but more generally ca. 25 K, whereas 20 have been studied employing 
pulsed Laval nozzles, which only generate temperatures as low as 53 K. Although there 
is certainly not a perfect overlap between the neutral-neutral reactions which have been 
studied to low temperatures using the CRESU method and those reactions whose rate 
coefficients are required for astrochemical modeling, the measurements have done much 
to identify the factors that control reaction rates at low temperatures, making it pos- 
sible to predict rate coefficients, or to extrapolate measured rate coefficients to much 
lower temperatures, with a fair degree of confidence. 

Based on a survey of the results for reactions between atoms or radicals with alkenes 
and alkynes, two criteria have been proposed (Smith et al. 20061 in order to decide 
whether reactions of this type will be fast; that is, have rate coefficients of ca. 10""'^'' 
cm'^ s~^ at 10 K. The first is if the observed rate coefficient at 298 K is > 10^^^ 
cm'^ s~^, the second is if the difference between the ionisation energy of the molecule 
and the electron affinity of the radical is less than 8.75 eV. These criteria have been 
supported by recent measurements on the reaction of 0("^P) atoms with a series of 
alkenes ( [Sabbah et al.|2007 l. 

Despite the CRESU and pulsed Laval studies, the experimental database of rate 
coefficients for reactions between neutral free radicals is sparse, especially at low tem- 
peratures. In general, the interaction between two radicals will give rise to more than 
one PES, as discussed above. Application of the correlation rules will show whether 
there are PESs that lead adiabatically to particular sets of products in specified elec- 
tronic states. In general, the two radical reactants may form a transient bond (an 
electron from each radical 'pairing up') and this interaction will lead to the lowest 
PES or possibly PESs. The lowest PES will probably exhibit a deep 'well' that can be 
accessed by a reaction path that exhibits no barrier. However, states of other multiplic- 
ities may be created and may not lead to products. One effect is that, if two radicals 
can undergo an exothermic reaction to yield two products, then the reaction is likely 
to occur rapidly, even at low temperatures. However, in estimating the rate coefficient, 
one should allow for the fact that some significant fraction of collisions will occur on 
PESs that either don't correlate adiabatically with the products or which do so over 
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PESs with significant energy barriers. This feature of radical-radical reactions necessi- 
tates the introduction of a factor < 1 which depends on the electronic degeneracies of 
the reactants and of the surface or surfaces over which reaction is judged to occur. 

In order to indicate how the values of rate coefficients for neutral-neutral reactions, 
and the uncertainties in these rate coefRcients, are chosen, we take as examples four 
"key" reactions listed later in the paper: those between (i) CN and NH3, (ii) CN and 
HC3N, (iii) O and C2, and (iv) N C3. 

The reaction between CN radicals and NH3 molecules was one of the first two 



neutral-neutral reactions examined at low temperatures by the CRESU technique ( Sims 



et al. 19941. The rate coefficient exhibits an unusually strong negative temperature- 



dependence; that is, it increases by a factor of more than 10 as the temperature is 
lowered from 295 to 25 K, the lowest temperature achieved in the experiments. The 
experimental value of the rate coefhcient was assigned a statistical error (at the level of 
95% confidence limits) of ca. 5%. As with most direct kinetic experiments on elemen- 
tary reactions, no information was obtained about the products of this reaction; the 
progress of reaction was followed by observing the decrease in concentration of the CN 
radicals. A recent theoretical analysis ( [Talbi and Smith|2009[ ), using a two-transition- 
state model and quantum chemical calculations of the minimum energy path for the 
reaction, reproduced the absolute value of the rate coefficient and its temperature- 
dependence, and concluded that the reaction proceeds exclusively to HCN and NH2, 
with zero yield of NCNH2 -|- H. This conclusion was subsequently supported by exper- 
imental measurements ( [Blitz et al.[2009) that demonstrate that the branching ratio to 
the latter set of products is certainly less than 5%. Almost uniquely, the rate coefficient 
for the reaction CN + NH3 NH2 + HCN at 10 K, fe = 5.0 x 10"^° cm^ s"\ can 
be assigned an uncertainty as low as A\ogiQ{k) = 0.15; this uncertainty arises mainly 
from the extrapolation of the experimental data from 25 to 10 K. The reaction to form 
NCNH2 + H is taken not to occur. 

The reaction between CN radicals and HC3N molecules represents the class of 
reactions where there are no low temperature measurements but there is an experi- 
mental determination of the rate coefficient at 298 K, in this case yielding fc(298 K) — 
1.7 X 10^ cm s~ (Halpern et al. ] |l989| . On the basis of both this rate coefficient 
and the relative values of the ionization energy of HC3N and the electron affinity of 
CN, the arguments of Smith et al. (20061 suggest that this reaction has no barrier 
and should be fast at low temperatures. The absence of any barrier to reaction is sup- 
ported by ab initio calculations ( Petrie and Osamura|20d4 |. The reaction is assigned a 
rate coefficient of 3.5 x 10"^*^ cm'^ s~^ at 10 K with an uncertainty corresponding to 
Zilog]^g(fc) = 0.6 at this temperature. 

The reactions of O atoms with C2 and N atoms with C3 are examples of reactions 
for which there are no experimental data in the literature. In the first of these reactions, 
CO and C atoms can be produced by a spin-allowed, highly exothermic reaction: 0(^P 
+ C2(^S+) CO(^S+) + C(3p); AH^gs = -480.8 kj mo\-\ By the criteria of 



Smith 



et al. ( 2006 1, this reaction would be expected to have no barrier. It is assigned values of 
k = 3.0 X 10"^*^ cm^ at 10 K and AlogiQ{k) = 0.5 at this temperature. The reaction 
between N atoms and C3 molecules in their (^Sj) ground state is problematical. 
Reaction to produce CN(^E^) + C2(^S^) is exothermic but spin-forbidden. However, 
the reaction producing C2 in its low-lying (^H) excited state is also exothermic but 
spin-allowed. On the other hand, N atoms are generally unreactive to unsaturated 
species, such as C2H2 and C2H4, in singlet electronic states, which is consistent with 
the rules proposed by Smith et al. (20061, because N atoms actually have a negative 
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electron affinity. Taking these considerations into account, this reaction is assigned 
a smaU rate coefficient at 10 K, fc = 1.0 x 10^^'^ cm'^ s^^, with large uncertainty, 
Z\logio(fe) = l-O- 



3.3 Radiative association reactions 



In association reactions, two atomic or molecular species, say A and B, combine to 
create one new molecular product, AB. In the ISM, such processes are important, 
since they allow the growth of species from smaller fragments. However, even in dense 
clouds, the gas density is too low to allow the coUisional stabilization of any initially 
formed adduct, which will have sufficient internal energy to re-dissociate, by collision 
with some other species. Consequently, the only mechanism by which the adduct can 
lose energy and generate stable AB is for it to radiate. The importance of radiative 
association has long been recognized: partly in respect of the reactions of C"*" and 



CHJ with H2 (|Herbst et al. 



1977 



Gerlich and Horning 19921. In both these cases 



the reactions forming two products (i.e., CH"'' + H and CHJ + H) are endothermic, 
so radiative association reactions are the only way to generate larger ions containing 
carbon and hydrogen. 

Although Table [2] lists separately the numbers of radiative association processes 
in the OSU model that involve (i) ion + neutral reactants, and (ii) neutral + neutral 
reactants, the mechanism for radiative association is the same whether or not one of 
the reactants is charged. The mechanism involves the initial formation of an energised 
complex by the association of the two reactants, followed either by the re-dissociation 
of the complex or the loss of sufficient energy by the spontaneous emission of radiation 
so that the complex then no longer contains sufficient internal energy to re-dissociate; 
that is 

A(+) +B ^ {AB(+'}* ^ AB(+) (14) 

In this scheme, the superscript (+) indicates that the reactant A may be neutral or 
charged. The superscript * indicates that the initially formed complex will possess 
sufficient internal energy for {AB*-^'}* to re-dissociate to the reactants A^"*"^ -I- B. 
A simple steady-state treatment of this mechanism shows that the second-order rate 
coefficient for radiative association, fcpjA, is given by: 



'^RA — '^ass { fcrad / C^rad + '^diss)} 



(15) 



In this expression, fcass is the second-order rate coefficient for the initial formation 
of {AbW}*, 

^rad is the first-order rate coefficient for the radiative stabilisation of 
{ab(+)}* , and fediss is the first-order rate coefficient for re-dissociation of the energised 
complex {AB*-"^^}* to reactants. This last rate coefficient depends on the strength of 
A(+) B bond and on the size of the system, that is the number of atoms in the {AB^'''''} 
complex. The more atoms, the more vibrations in {AB^^-*} in which the energy can 
be distributed, and the slower re-dissociation will be. fcjisg also depends strongly on 
the internal energy in {ab(+)}* or, in thermal systems, the temperature, so that rate 
coefficients for radiative association are expected to show a strong negative dependence 



on temperature. The equation (151 is for a thermal (canonical) system; more detailed 
microcanonical formulations using phase space theory are also available ( [Bates and| 
19881. For systems in which there is a competitive two-body exit channel, 



Herbst 



the denominator of eq. ( 15 1 must also include a term for adduct dissociation into 
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products. This term need not reduce the rate coefficient of radiative association to an 
immeasurably small value if there is a barrier in the exit channel that lies below the 
energy of reactants. It is also possible that association and normal exothermic channels 
can occur via different, unrelated reaction pathways. 

For small species, k^^gg » ftj-^d, and radiative association is an inefficient process; 
that is, it occurs in only a small fraction of the collisions in which {AB^^^}* complexes 



are formed. In these circumstances, eq. ( 15 1 simplifies to 



'i^RA = fcrad{fcass/fcdiss} (16) 

Canonical treatments using this equation can make use of the approximation that 



fcass/fcdiss is a thermal equilibrium constant (Bates and Herbst 19881. Moreover, in 
interstellar environments, radiative association is only a significant process in reactions 
where either one reactant is the (relatively) abundant H2 or when {AB*-^'}* is 'large': 
that is, it has more than ca. 10 atoms, so that k^^ss < fefad and fcj^A ~ ^ass. 

The most difficult part of estimating rate coefficients for radiative association - 
and this applies to ion-neutral as well as neutral-neutral systems - is generally es- 
timating fcrad- In those cases where only the ground electronic state of the 
species is involved, the radiative process comprises transitions from vibrational levels 
in {AB^"*"^}* to vibrational levels in AB''~''\The rate of radiative emission in this case 



can be obtained from the equation ( |Herbst|1982 1 

1 

ME) 



modes 

1 . 



ferad = TT^ 2^ ^',1,0 2^ niPs-i{E - nihvi) (17) 



= 1 ni = l 



where E is the internal energy in the complex; ps{E) and ps-i{E -— nihvi) are the 
density of states associated with s oscillators at energy E and with s - 1 oscillators 
at energy {E — riihvi), the energy nihvi being 'locked' in the ith oscillator. Values of 
Ai l 0: Einstein coefficient for spontaneous emission from the = 1 to = level 
of the ith oscillator, can frequently be obtained from theory or, in some cases, from 
experimentally determined infrared absorption coefficients. Values of fc^ad are generally 
found to lie within the range 10 - 1000 . Because the spontaneous emission coeffi- 
cients depend on the cube of the emission frequency, higher values of k,.ad are found 
for molecules containing hydrogen atoms because they will possess some relatively high 
vibrational frequencies. One major source of uncertainty in estimating rate coefficients 
for radiative association arises from the possible role of electronically excited states of 
If such states correlate with 

A(+) 

+ B, then the value of fc^ad may be dramat- 
ically increased, even if electronic transitions between the electronically excited state 
and the ground state are not fully allowed. Such a possibility occurs for the C+ + H2 
system ( [Herbst et al.|1977 l. 



Although similar mechanisms for radiative association operate, whether or not one 
of the reactants is electrically charged, laboratory measurements - indeed laboratory 
measurements at very low temperatures - have been possible only for ion-neutral sys- 
tems, but not for neutral-neutral reactions. Ion-neutral radiative association can be 



studied down to very low temperature using ion traps (Gerlich 20081. The reaction 
C*" -I- H2 — >■ CHJ -I- hu has been studied using both n-H2 and P-H2. At 10 K, the 
rate coefficient is 7 x 10~^^ cm^ s~^ for n-H2 and 1.7 x 10~^^ cm^ s~^ for P-H2, with 
uncertainties of only ca. 50%. These results are in reasonable agreement with theory 
( Gerlich and Horning|l"992 1. 
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Especially in the case of the radiative association between neutral species, where 
there is little in the way of direct experiments, the rate coefficients have to be estimated 
in one of two ways, both of which require a theoretical estimate of the rate coefficient, 
fci-adi for the radiative emission step. First, one can make use of experimental data on 
the related process of coUisionally-assisted association. In these reactions, the radiative 
stabilisation step in the mechanism represented by eq. (15) is replaced by one in which 
energy is removed from the energised complexes, {AB'-"*'-'}*, in collisions with a 'third- 
body species (M)'. Then eq. (16) for k^A can be replaced by 



fe3rd[M] = fcM[M]{A:ass/fcdisJ (18) 



Here kyi is a second-order rate coefficient for the coUisional process in which energy 
is removed from {AB*-^-*}* and k^,-^ is a third-order rate coefficient for coUisionally- 
assisted association under conditions where k^[ss » ^m[M^]- It should be mentioned 
that the implicit assumption is made that individual collisions can stabilize the adduct; 
in reality, master equation techniques, in which many collisions might be needed, are 
preferable. Secondly, when there is no experimental data on the coUisionally-assisted 
association process, the ratio {fcass/fcdissli which corresponds to an equilibrium pop- 
ulation of internal states in the energized complex, can be estimated using thermal, 
quasi-thermal, and phase-space approaches ( [Bates and Herbst|1988[ ) as long as there is 
structural and vibrational information about the adduct. In general, theoretical treat- 
ments are at best accurate to an order-of-magnitude, for both ion-neutral and neutral- 
neutral systems. If there is a larger disagreement for a system that has been studied 
experimentally, it is likely that the structure of the theoretical adduct is incorrect or 
that the adduct is stabilized electronically. Use of three-body data can also pose an 
additional problem, since many experiments are performed in a near-saturated limit, 
in which the third-order rate coefficient obtained is only a lower limit to the true value. 

An association process in which a neutral species with a large electron affinity, 
typically a radical, and an electron stick together to form an anion via emission of a 
photon is known as radiative attachment. This reaction is thought to be the primary 
one for the production of the interstellar anions - CaN", C4H~, CgH", and CgH" - 
observed to date ( Herbst and Osamura|2008 1. At present, there have been no measure- 



ments and no detailed calculations of radiative attachment rate coefficients for these 
systems. A statistical theory to determine the rate coefficient, akin to that for radiative 
association, was developed by ( Herbst|1981 \ and can be expressed in a form similar to 



the radiative association rate coefficient in eq. (|15f . Here the attachment proceeds by 
way of a temporary adduct, which radiates away energy as in radiative association (see 



eq. ( 17 1). If it is assumed that the formation of the adduct occurs with s-wave electrons 
only, phase-space theory can be used to derive a simple form for fcatt. Use of this theory 
shows that the larger detected interstellar anions are formed on nearly every collision 



between radicals and electrons, with a rate coefficient in excess of 10 ^ cm'^ s ^ (iHerbst 



and Osamura|[2008 1 . There is, however, a large (two order-of-magnitude) discrepancy 



between model and observations when the theory is used for the radiative attachment 



of C4H (Harada and Herbst 2008, ). More detailed calculations are needed. These may 



follow the very detaile d treatment of coUis ional and radiative electron attachment of 
SFe which forms SFg \Troe et al.|2007a|b I. 
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3.4 Dissociative recombination reactions 

Dissociative recombination (DR) denotes the process whereby a molecular ion recom- 
bines with a free electron and the energy released is dispersed through dissociation of 
the intermediate neutral into two or more fragments, with conversion of any excess 
energy into kinetic and internal energy of the fragments. The importance of DR is 
two-fold: firstly, it can serve as the final step in the synthesis of neutral molecules by 
transforming a protonated species, which was previously produced in a series of ion- 
neutral reactions, into the unprotonated species. Such reaction sequences have been 
invoked as synthetic pathways for a multitude of different interstellar molecules; e. g., 
HCl plake et aLpOSS] ): 

H2 + C1+ ^ HC1+ + H, (19) 
H2 + HC1+ H2CI+ + H, (20) 

H2CI+ + e" ^ HCl + H. (21) 

In this mechanism, reaction (21) corresponds to the DR. Secondly, protonation followed 
by DR can serve as a very efficient destruction mechanism for molecules with high 
proton affinities, such as nitriles. Several ions themselves, such as cyclic CaH^, are 
sufficiently unreactive in their environments that DR constitutes their sole destruction 
process ( Keller et al.||1998 l. DR reactions therefore form an integral part in chemical 



networks of interstellar clouds . Since the abundances of important species predicted 
by models can be greatly affected by the rates and product branching fractions of DR 
reactions, knowledge of these parameters is vital for accurate predictions. 

Thermal rate coefficients of DR processes vary, but not to an excessively large ex- 
tent. In fact, most of them lie in the range between 1 x lO^'' and 4 x 10~® cm'^ s~^ 
at 300 K, with larger ions tending to show higher rate coefficients. However, many ex- 
perimental measurements of DR rate coefficients were performed at room temperature 
and it is not entirely trivial to extrapolate these parameters to the low temperatures 

encountered in cold dense cores. As a barrier-less electron-ion reaction, DR reactions 

1/2 

should display a temperature dependence of the thermal rate coefficients of T ' , 
which is also assumed in models if no contrary evidence exists. Nevertheless, as experi- 
ments have shown, deviations from this behavior are common: usually the temperature 
dependence is stronger, leading to higher values of k at very low temperatures for a 
variety of reasons including resonances and potential crossings. 

Branching fractions of DR reactions have proven to be surprising and often coun- 
terintuitive to chemical common sense. According to an early theory formulated by 



Bates (19861, DR reactions should preferably follow the product pathway involving 
the least rearrangement of valence orbitals. In the case of protonated neutrals, this 
implies formation of the unprotonated parent substance and a hydrogen atom. Al- 
though so me ions such as HC O"*" follow that rule, many other important species such 
as H3O''' ( Jensen et al.|[2000 l and protonated methanol (CHaOH^) do not ( ^Geppert 



et al.|2006| . In the DR of these species, break-up processes of the intermediate neutral 
into three fragments dominate, which should be highly unlikely according to Bates's 
theory. Another early theory, that of |Herbs"t] ( |1978[ ), uses phase-space theory and can 
only handle three-body channels if they occur via breakup of a primary species formed 
in a two-body channel. This theory is reasonably successful for the products of HsO^ 
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+ e~ but is in general inaccurate. Therefore one has to rely on experimental measure- 
ments, which are hitherto unavailable for many ions, so many branching ratios used in 
models are still educated guesses. 

Several experimental methods have been employed to determine the rate coeffi- 
cients and branching ratios of DR reactions. The first studies were carried out using a 
stationary afterglow (where the ions are produced in a discharge followed by subsequent 
ion-neutral reactions to form the desired ionic species, thus in the " afterglow" of the dis- 
charge) ( Biondi and Brown|1949 1 . The faults in this approach were partially overcome 
by the more novel flowing afterglow method. Flowing afterglow methods have been 
employed with success for many systems. However, they also have their limitations. 
Most studies are carried out at room temperature and therefore, as discussed above, 
have limited validity for the interstellar medium, although recently experiments on the 
DR of H!^ have been performed at 77 K ( Glosfk et al. |2009| . Furthermore, it is difficult 
to determine the branching ratios of a DR reaction, since rarely can all the product 
channels be quantified. These facts led to the employment of merged beam techniques 
to investigate DR. Magnetic storage rings, like ASTRID (Aarhus, Denmark), TSR 
(Heidelberg, Germany) and CRYRING (Stockholm, Sweden), have proven especially 
successful. However, experiments in storage rings do have their shortcomings. Firstly, 
there is usually no real control of the rovibrational states of the ions, although in the 
case of small ions (e. g. H^) calibrated ion sources have been proven to produce ions 
with a distribution between different states corresponding to temperatures relevant to 
the ISM (McCall et al. 20031. Secondly, mass selection does not allow separation of 



isomer-pairs such as HCO+ and HOC+, which can be produced simultaneously in the 
ion source and could show different behavior upon DR. The fact that both storage ring 
and afterglow methods have their merits and faults is painfully illustrated by the fact 
that for the DR of CH^ ion storage ring measurements yielded a predominance of the 



CH3 -I- 2H channel ( Zhaunerchyk et al.|2009 l, whereas flowing afterglow studies found 
the CH4 -I- H channel to dominate ( Adams et al.|2009J . Imaging measurements showed 
that the CH3 + 2H production proceeds in a two step mechanism, flrst leading to H 
and a highly excited CH4 molecule, which could subsequently dissociate to yield CH3 
-|- H. At the higher pressures in the flowing afterglow experiments, the CH4 interme- 
diate can be deactivated by collisions with buffer gas, leading to a higher branching 
fraction of the CH4 -I- H pathway. 

Given the problems associated with experimental determinations of the branching 
ratios and rate coefBcients of DR reactions, one could wish to resort to predictions by 
high-quality ab initio calculations. However, due to the involvement of highly excited 
states and complicated potential surfaces, such computations are very challenging. 
Also, a large number of open electronic states exist, which might be close in energy, 
bringing about break-down of the Born-Oppenheimer approximation. Furthermore, not 



only bound states, but also free electron states have to be included (Larsson and Orel 



2008 1, making such studies very time-con suming. Although the oretical investigations 
have been performed for the DR of HCO"'' ( Mikhailov et al.|2006 1, it is quite reasonable 
to assume that such efforts will , at least in the foreseea ble future, be limited to rather 
small systems, such as for Hl^ (Dos Santos et al. |2007| . 

Finally, it should be mentioned that, if electrons undergo radiative attachment 
efflciently to large neutrals such as PAH's, it might be possible that the abundance 
of negatively-charged PAH's exceeds the electron abundance. In this case, it would be 
necessary to include mutual neutralization processes between small cations and large 
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Fig. 1 Flow diagram of Uncertainty Propagation and Sensitivity Analysis 

anions in interstellar models. Unfortunately, the amount of information concerning 
such processes is miniscule. More laboratory studies are certainly needed. 



4 Uncertainty analysis, uncertainty propagation and sensitivity analysis 

Because chemical models contain hundreds to thousands of chemical reactions, the 
model predictions can be strongly affected by uncertainties in the rate coefficients. 
Dealing with uncertainties in modeling is a three-stage process: 

1. Uncertainty Analysis (UA) consists in identifying the uncertain inputs of the model 
and quantifying their uncertainty. The type of information to be issued at this stage 
is partly conditioned by the next stages: it can be a simple standard deviation or 
a full probability density function; 

2. Uncertainty Propagation (UP) is used to transfer the uncertainty in the model 
input(s) to the model output(s); 

3. Sensitivity Analysis (SA) enables one to identify input (key) parameter(s) with a 
large contribution to the uncertainty in the model output (s). 

The methodological framework presented in Fig. [l] according to |de Rocquigny et al.| 



(20081, summarizes the operations and information flow involved in a full analysis. 



4.1 Uncertainty Analysis 

UA is challenging for the modeller, because it often requires proficiency in the ex- 
perimental techniques used to evaluate the quantities of interest. A general problem is 
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that the reference databases generally lack quantified uncertainty statements about the 
stored data, or at least properly quantified for the purposes of UP/SA. As discussed 
in this article, KIDA is an attempt to address this shortcoming. 

SA results are quite sensitive to this elicitation step: too small uncertainties lead to 
an over estimation of prediction precision, whereas too large uncertainties can lead to 
useless predictions. At the present stage, we are still in need of clearly defined guiding 
rules, and groups of experts, containing both experimentalists and modellers, to work 



on this topic, as has been done in respect of the Earth's atmospheric chemistry ( Sander 

|et al.|2006| [Atkinson eraL][2004|. 

Following the work of Dobrijevic and Parisot ( 1998 1 and Dobrijevic et al. ( 2003 1 on 
giant planets, exhaustive analysis of these issues was started for Titan's atmosphere. 



An evaluation of the uncertainties for neutral-neutral reactions (Hebrard et al 



20071 and for ion-neutral reactions (Carrasco et al. 2007 Carrasco and Pernot 



2006 



2007 



Carrasco et al.||"2008 1 has been undertaken, in order to prioritize reactions in need of 



better characterization, possibly by laboratory measurements. Application to Titan's 
atmosphere shows that uncertainty in the chemical rate coefficients strongly reduces 
the predictivity of atmospheric models (Hebrard et al. 20071. A particular issue in 



this system (which is also relevant to astrochemistry) is the use of extrapolation for 
90% of the reaction rates to the low temperatures of Titan's atmosphere (50-200 K). 
The task of providing reasonable uncertainty on extrapolated data is nevertheless an 
unavoidable step in order to prioritize key reactions, in need of better determination 
fHebrard et al.|2009 l. 

A similar approach has been undertaken by Wakelam and coworkers for astrophys- 
ical objects | | Wakelam et al.||2005[ |2006| [Vasyunin et al.||2008| [Wakelam et al.||2009[ ), 
as presented in the following sections. 



4.2 Uncertainty propagation 

Uncertainty propagation is a tool to quantify the precision of model predictions. Be- 
cause of their central role in all fields of science and engineering, UP methods have 
been standardized at the international level ( [BIPM et al.[|2008a[ [BIPM et al.[|2008b l. 
They are of two kinds: (i) combination of variances, as described in the Guide to the 
Expression of Uncertainty in Measurements (called "the GUM") ( [BIPM et al.|2008a[ ); 
and (ii) propagation of distributions, as described in the first supplement to the GUM 
(GUM-Suppl) ( [BIPM et al.|2008b| . 



Denoting as the i input of a model G, and Z as an output of this model 
{Z = G{Xi, X2, ■ ■ ■ , Xn)), the GUM stipulates that uz, the standard uncertainty of 
Z, should be calculated by the following formula: 



2 



X 



dG 
dXi 



cov{Xi, Xj 



(22) 



where cov{Xi,Xj) is an element of the covariance matrix of input parameters, with 
variances as the diagonal elements coY^Xi, Xi) = u^x ■ The so-called sensitivity coeffi- 
cients {dG / dXi)^ are evaluated at the mean value of input parameters X. 

This method is exact if the model G is a linear function of the inputs, and if the error 
distribution of these inputs is symmetric. For small relative uncertainties (typically less 
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than 30%), it provides generally an excellent approximation of uz, even for non-linear 
models. 

However, astrochemistry and atmospheric chemistry have to deal with uncertain- 
ties much larger than 30% ( Hebrard et al.|2006 Vasyunin et al.|2008 l, and because of 
the non-linearity of photochemical models, eq. ( 22 I cannot be relied upon. It neverthe- 
less reveals that there are two essential components to uncertainty propagation: the 
sensitivity of the model to the inputs; and the amplitude of input uncertainty, which 
should be present in any other UP method. 

The method of propagation of distributions requires the definition of probability 
density functions (pdf) for all the parameters (BIPM et al. 2008b I. The reference 
equation for the propagation of distributions is 



p{Z^a)= / dXidX2...dX„5(a-G(Xi,X2,...,X„)) p(Xi,X2,...,X„ 



(23) 



where a is a value in the range of Z. This equation provides the pdf of Z , p{Z) 
by transforming the pdf of the inputs, p{Xi,X2, ...,X„), through the model. It can 
rarely be solved analytically. Instead, the GUM-Suppl recommends the use of Monte 
Carlo Uncertainty Propagation (MCUP) methods (see also Section [5|. The principle is 
simple: one generates a number A'run of random input sets from their pdf, and for each 
set, one evaluates the value of the output, thus generating a representative sample of 
the output pdf p{Z). This sample is then used to estimate quantities of interest. 

In the hypothesis of independent inputs with Normal/Gaussian type uncertainties, 
the same ingredients are necessary to generate the input random sample as for the 
method of combination of variances; i.e., average values and standard uncertainties. 
Other distributions might require more complex elicitation procedures, especially in 
the case of correlated inputs {e.g. branching ratios Carrasco and Pernot|2007[). 



4.3 Sensitivity analysis 

Sensitivity Analysis (SA) methods are not as standardized as UP methods: a variety 
of methods are more or less adapted to specific scenarios, but globally, SA consists in 
the apportionment of output variance between the uncertain inputs ( Saltelli et al.|2000 
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In this work, we are interested in the identification of key parameters being respon- 
sible for a large part of the output variance. In the case of a linear model, this could 
be done by comparing the contributions of inputs in eq. |22[ The ingredients are thus 
the model sensitivity factors and the input variances/covariances. Combination of both 
terms means that an output might be insensitive to a very uncertain input, or very 
sensitive to a fairly well-known parameter. This remains true for nonlinear models, and 
both contributions remain essential. Sensitivity analysis has therefore to be based on 
careful uncertainty analysis and uncertainty propagation steps. 

A convenient SA method can be combined with MCUP through the calculation of 
nput /output correlation coefficients ( Helton et al.|[2006 K Pearson's correlation coeffi- 



1 

cient 



C(Xj , Z) = ^ _ = (24) 

'Eiix',-x,rj:^{z^-zr 
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is often used ( [Wakelam et al.|2005| |2006| |Vasyunin et al.||2008| [Wakelam et al.||2009 1 



where the sums run over the sample length (1 = 1, Nmn)- Correlation coefficients vary 
between -1 and 1, and influential parameters are identified through the comparatively 
large absolute values of their correlation coefficient with the output \C{Xi, Z)\. Spear- 
man's rank correlation coefficient is considered more appropriate for nonlinear models 
( |Hehon et al.|2006| [Dobrijevic et al.|2010| ) , albeit our experience with sensitivity anal- 
ysis of chemical networks is that the differences between both methods are generally 
minor. 

Other SA methods, based on specific sampling patterns, can be more economical to 
run if one does not require UP results ( Saltelli et al.|[2000 l. In the present work, both 



MCUP and SA stages are of interest, and the Correlation Coefficients SA method is a 
computationally inexpensive byproduct of MCUP. 

In the context of sensitivity analysis, a key reaction has a reaction rate uncertain 



enough to contribute significantly to the model output uncertainty ( 


Wakelam et al. 


2005IJ2006 


IVasyunin et al.|2008| |Dobrijevic et al.|2010||Hebrard et al.|2009| Wakelam 


et al.|2009 


. Identification of such key reactions enables one to point out those in need 



of further study, either by an evaluation committee, or by new experiments or calcula- 
tions. The improvement of the precision of model predictions is a slow iterative process 
between stages of sensitivity analysis and database updates, as new key reactions are 
expected to emerge when the previous ones have been updated ( Hebrard et al.|2009 l 



In order to select key reactions, a target has to be defined. This could be the 
abundance of a single species of interest. In that case, one evaluates and compares the 



input /output correlation coefficients for a single output species (Wakelam et al.||2005 



Dobrijevic et al.|2008[ ). If one is interested in a global improvement of model precision 



an additional criterion is necessary, as for instance selecting reactions presenting large 
correlation coefficients with a large number of output concentrations (I Wakelam et al.l 



2005 2006 Vasyunin et al.||2008 Hebrard et al. 112009 Wakelam et al.||2009). The key 



reactions can then be ranked according to the number of species they affect. 



5 Gas-phase molecular cloud chemistry 

To improve the chemical modeling of cold cores, we applied sensitivity methods to 
identify " key reactions" . A detailed review of the existing experimental and theoretical 
data on these reaction rate coefficients was then conducted to propose new values or 
confirm the ones already used in the models. We also studied the impact of the proposed 
new rate coefficients on the cold core chemistry. These three steps are described below. 



5.1 Sensitivity analysis: identification of key reactions 

To identify key reactions, we used Monte Carlo sensitivity methods as described in 



Wakelam et al. ( 2009 1 and this review. Briefiy, we randomly varied all the reaction 



rate coefficients within a certain factor. For each set of rate coefficients, we computed 
the chemical evolution of the cloud and thus obtained A'' chemical evolutions. For each 
species j and each chemical reaction i, we calculated Pearson correlation coefficients 
(see Section [431 Eg. [24]): 



P;(f) = C(log(A,(t)),log(fe,)), 



(25) 



24 



-6 




log(time [yr]) 

Fig. 2 Abundance of HC3N as a function of time. Each black line is the result of one run 
(one set of reaction rate coefficients). The density is 2 X 10* cm~^, species are initially in the 
atomic form (except for H2) and the rate coefficients are varied within their uncertainty range, 
as defined in the database. 



where Xj{t) is the abundance for species j and ki is the rate constant of reaction i. 
The larger the value of PJ, the stronger the correlation and the more important the 
reaction i for the abundance of the species j. 

For this work, we used the nahoon chemical modelj^with the osu-09-2008 network. 
In addition, we allowed gas-phase species to stick onto dust grains and to be evaporated 
from the surfaces by thermal desorption and cosmic-ray heating. We assumed a fixed 
homogeneous temperature of 10 K, a visual extinction Ay of 10, and a cosmic-ray 
ionization rate (^h2 of 1-3 ^ 10~^^ s~^. We considered two densities: 2 x 10** and 
2 X 10^ cm^"^, which refer to the total nuclear density of hydrogen. The typical low- 
metal elemental abundances were used (see Table 1 in |Wakelam and Herbst|2008| ). We 



considered two types of initial conditions. In the first, all the species are in the atomic 
form except for hydrogen, which is entirely molecular, while in the second, half of the 
ionic carbon is in CO, and the other half is still in C"*". 

To sample the reaction rate coefficients, we used two different sets of uncertainty 
factors. In the first set, the uncertainty factor is defined in the database, where most 
reactions are assigned a factor of 2, as originally suggested in the UDfA database, 
but others, studied experimentally, have a smaller factor of 1.25 or 1.5 All radiative 
association reactions have a larger uncertainty factor of 10. In the second set, the 



http ; //www.obs.u — bordeauxl.fr/amor/VWakelam/ 
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uncertainty factor is fixed at 10 for all reactions. The goal of this second set is to see 
the effect of varying the rate coefficients over a larger range of values. 

Considering the different densities, initial conditions and uncertainty factors, we 
constructed 8 different models and for each of these models, we made 2500 runs by 
randomly varying the rate coefficients using a log-normal distribution. As an example, 
we show the results for HC3N in one model (an H density of 2 x 10* cm~"^, atomic 
initial conditions, and the uncertainty factors as defined in the database) in Fig. |2] 

In Table [Sj we report the reactions that have a significant impact on more than 
one species. To determine these reactions, we count for each reaction the number of 
species with a correlation coefficient larger than 0.3 and an abundance larger than 
10""'^^. The abundance criterion is chosen because these species can be observed in 
dense cold cores and thus used to put constraints on the model. We used two different 
times: 10^ and 10^ yr, covering the probable age of most cold cores. A more thorough 
approach would be to use a larger number of times, because some reactions important 
at earlier stages may contribute to the synthesis of the larger molecules at the time 
chosen. We use a Boolean approach: either a species contributes to the importance of 
a reaction or not. Another possibility, which is considered in Section 6, is to sum the 
correlation coefficients for all species. Then, the larger the sum, the more important 
the reaction. In both cases, the result depends on the chosen uncertainty factor in the 
rate coefficient. Each method has its pros and cons. In this current study, we checked 
that the result of the analysis is not significantly changed if the second approach is 
used. The list of key reactions is very similar for the different models; only the relative 
importance of the reactions really changes. 

Tables [3] and [4] give the key reactions in the following case: density of 2 x 10** cm"'^, 
species initially in the atomic form (except for H2) and rate coefficients varied within 
their uncertainty range, as defined in the database. The 7 other models we ran gave 
the same reactions but in a different order of importance except for the first reaction in 
TablejS] The association reaction C + H2 — )■ CH2 + hiy, which heads the list, influences 
about 73 abundant species at 10^ yr. Although association reactions can have small 
rate coefficients (see Section [3.3[ ), if the reactants are abundant, these reactions can 
play an important role. In Table |4j we show the key reactions, divided as to type, for 
individual observable species. Here again, we considered only species with a correlation 
coefficient larger than 0.3 and an abundance larger than 10~^^ at 10^ or 10^ yr. 

5.2 Reaction review 

Of the 59 reactions that have been identified as particularly crucial in the chemistry 
of cold cores in Tables |3] and |4j 34 have been reviewed and new values of the rate 
coefficients are proposed for 22 of these reactions. The new rate coefficients that we 
propose as a result of this review are compared in Tables [5] [6] [7] and [8j according to 
type of reaction, with the values that are currently in the OSU and UDfA databases. 
We have considered rate coefficients for temperatures between 10 and 298 K. In the 
tables, a designation of zero means that the rate coefficient is immeasurably small. 

The results of this effort have been included in a new Kinetic Data base for Astro- 
chemistry (KIDA), which is available at http://kida.obs.u-bordeauxl.fr'. This database 
differs from the OSU and UDfA data bases in that, at least for the key reactions 
identified in this article, the reasons for the choice of particular values of the rate co- 
efficients and their uncertainties are explained in Datasheets (accessed by clicking on 
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Table 3 List of key reactions by number of species influenced. 



Reaction 



Strongly affected species^ 



CH+ 



H2 CH2 + hi/ 



C2H+ 



H2 ^ 
f H2 
CH+ + CO 
C2H+ + e- 
HSiO+ 
HSiO+ 
C3H+ 



C3H+ 
CH+ 4 
C4H+ 
CH+ 4 
C4H+ 



■f e- 
f e- 

- H2 - 

- H2 - 
HCN 

- H ^ 



CHJ + hu 

C2H+ + hu 
► C2H3O+ + hiy 
* C2H3 + H 
-s> SiO + H 
-s> Si + OH 
^ C3H+ + hu 
■> H3CI + hu 
C2H4N+ + hu 



C4HJ + hu 
NH3 -5- CH6N+ + hu 
f O -s> HC4O+ + H 



73 
18 

C2H2O, C2H3, C2H+, C2HO+, C2H2N+, C2H+ 
C2H2O, C2H3O+ 
C2H2O, C2H3 
Si, SiO 
Si, SiO 

C3H+, C3H+, C3H+, H3C+ 
C3H+, C3H^, C3H+, H3C+ 
C2H2N, HC3N, C2H3N, C2H4N+ 
C4H2, C5H, CeHe, C4H+ 
CH3N, CH5N 
C3O, HC4O+ 



C3H2 
C3H2 



H2C3, 
H2C3, 



^ Right column indicates the number or names of the species influenced strongly by each reac- 
tion. 




Chemical models 

(dense clouds, disk, Titan...) 



T 



Sensitivity 
Methods 



Group of experts in 
chemistry 



Review and propositions 
^ for reference values ^ 



Important reactions 
(depends on the system and the 
chemical network) 



Fig. 3 The construction of the KIDA database. 



the entry 'Details' in the first page devoted to a particular reaction). An example of 
datasheet for the reaction iT^ + O is given in Appendix |a| Figure [s] illustrates 
the construction of this database. 

These DataSheets follow the form adopted by the lUPAC Subcommittee for Gas 
Kinetic Data Evaluation, which reviews and evaluates kinetic, photochemical and het- 
erogeneous data for atmospheric chemistry ( http:/ /www. iupac-kinetic.ch.cam.ac.uk) . 
For each reaction the Datasheet includes: 

1. the possible product channels, where appropriate, and the electronic states of the 
reactants and products; 

2. the enthalpies of reaction for the possible reaction channels; 
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Table 4 List of key recLctions for specific species 



Neutral-Neutral reactions Affected species 

C + C3O ^ Cs + CO cio 

C + OCN CO + CN OCN 

H + CH2 CH + H2 CH 

O + CN -5> CO + H CN 

N + CN -5> C + N2 CN 

O + NH -5> NO + H NH 

O + C2 -5> CO + C C2 

O + C2H -5> CO + CH C2H 

O + C3H ^ C2H + CO C3H 

N + C3 ^ CN + C2 Ca 

N + NO N2 + O NO 

O + NH2 ^ HNO + H NH2 

O + HNO NO2 + H N2O 

O + HNO N2O + H N2O 

CN + NH3 -> NH2CN + H NH2CN 

O + C3N -5> CO + C2N C3N 

N + C4N ^ CN + C3N C4N 

N + C4H C4N + H C4N 

N + C2N ^> CN + CN C2N 

CN + HC5N NCeN + H NCeN 

CN + HC3N NC4N + H NC4N 

Association reactions Affected species 

C4H+ + HC3N -> C7H3N+ + hi/ HC7N 

HS+ + H2 ^ H3S+ + hi/ H2S 

HCO+ + H2O CH3O+ + hi/ CH2O2 

C+ + H2 ^> CH+ + hi/ C3 

S + CO -5> OCS + hu OCS 

CH+ + HC3N -> C4H4N+ + hi/ CH3C3N 

CH+ + HC5N ^ C6H4N+ + hv CH3C5N 

Si+ + H2 ^ SiH+ + hi/ HNSi 

Ion-neutral reactions Affected species 

C2H+ H- O C2H2O+ + H C2H2O+ 

C+ + S ^ S+ -I- C H2CS 

CbH+ -f N C6N+ -1- H C5H2N+, HCbN 

C2H+ + N ^ C2H2N+ + H2 C2H2N+ 

C4H+ + S ^ HC4S+ -I- H C4S 

H+ C ^ CH+ + H2 C3H3 

HJ -I- O ^> OH+ + H2 H3O+ 

C2H+ -I- N ^ C2NH+ -I- H2 C2NH+ 

Dissociative recombination Affected species 

HC4O+ -I- e- ^ C3O -I- CH C^O 

H2NC+ + e- ^ HNC + H HNC, H2NC+ 

C5H2N+ + e- C5N + H2 C5H2N+ 

HC4S+ + e- C4S -f H C4S 

H2CO+ -h e- CO -h H -t- H H2CO+ 

CNC+ -h e- ^ CN -h C CNC+ 

HC4S+ -t- e- ^ C4S -t- H C4S 

HC3S+ -f e- -)■ C3S H- H C3S 

HC3S+ +e- C2S + CH C3S 
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3. a Table giving experimental (and theoretical) values of the rate coefficients, and 
the values given in previous reviews, including the OSU and UDfA data bases; 

4. comments on the experimental and theoretical studies that have been performed; 

5. recommended values of both the rate coefficients and their estimated uncertainties 
for temperatures between fO and 298 K except as mentioned in the text, and 
comments (and justification) for the choices that have been made. 

When we have an estimation for the uncertainty factor Func, we have indicated it in 
the Tables [5] [6] [7j and|8]otherwise, the symbol - is used. There are three reactions that 
have been examined - C + H2 — >■ CH2 + hi/ and C"*" + S — >■ + C - but for which 
we do not give recommendations at the present time, since no experiment or detailed 
calculation has been done at low temperature. The impact of their rate coefficients on 
the abundances of the species appears to be especially strong. They will be discussed 
separately in the next section. Finally that the association reaction S + CO — )■ OCS + 
hv is currently being studied theoretically by the group of Gunnar Nyman (Goteborg 
University, Sweden). 
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2 3 4 5 6 7 8 

log(time[yr]) 

Fig. 4 The D norm parameter (see text) as a function of time when the different types of 
reactions are updated: DR stands for dissociative recombination; RA, for radiative association; 
Neutral, for neutral-neutral reactions; Ion, for ion-neutral reactions; and All, for the case in 
which all rate coefficients are updated. 

5.3 Impact of new rate coefficients 

We separate this discussion into four parts. In the first part, we show the impact of the 
updated reactions listed in Tables[5]|6][7|and[8] We compare the computed abundances 
with observations in dense clouds in the next part. Then two reactions, C -I- H2 — > CH2 
-I- hu and + S — !> C + , will be considered separately because a correct evaluation 
of their rate coefficients at low temperature requires detailed calculations which have 
not been carried out as part of the present study. We will however demonstrate how 
sensitive the abundances of the species are to the rate coefficients of these reactions 
assuming reasonable minimum and maximum values. For this section, we will use the 
following physical conditions: a temperature of 10 K, an H density of 2 x 10** cm"'^, 
and an Ay of 10. 

5.3.1 Impact of updated reaction coefficients 

In order to study the impact of the new rate coefficients, we have calculated at each 
time a function Dnorm{t), defined by the equation 

D^or^it) = 1 ^ \log{Xi,^{t)) - log(X^pd,,,(t))l (26) 
j 
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Fig. 5 O and H2O abundances as a function of time predicted by the model in the case of the 
osu-09-2008 network (solid line) and when the ion-neutral reaction + O rate coefficient 
has been updated (dashed line). 



with X^^^ the abundance of species j computed with osu-09-2008 and -'^updato^^) 
abundance of species j computed with the updated networks. The sum A'^, which is 
102 for aU networks, is over aU the species whose abundance is larger than 10""'^^ at 
10'^ yr. In addition to the osu-09-2008 network, we have used 5 other networks. In four 
of the five networks, we have updated separately each class of reaction: neutral-neutral 
reactions from Table |5j radiative associations from Table [6] ion-neutral reactions from 
Table [7| and dissociative recombinations from Table [8] In the fifth network, we have 
updated all the types of reactions at the same time. 

The parameter Dnorm is shown in Fig. [4] for the different networks. The changes 
in the rate coefficients for the radiative association reactions are more important for 
times earlier than 10^ yr whereas ion-neutral and neutral-neutral reactions impact the 
abundances of the species especially between 10^ and 10^ yr. In particular, there is a 
peak of Dnorm at 5 X 10^ yr due to the change in the rate coefficient for the -|- 
O reaction (since it is the only updated ion-neutral reaction) that affects strongly the 
abundance of many species including C, C''', O, O2, SO, CH2 and II2O at that time 
(see Fig. [5|. This reaction is indeed the main destruction path of the atomic oxygen. 
The abundance of O is then decreased whereas the ones of H2O and O2 are increased. 
The update of the branching ratios does not affect strongly the results at steady-state 
since OH"'' is rapidly converted into 1120"*" by reacting with II2. 

When neutral-neutral reactions are updated, the species that are most affected 
are C2, NO, C2H, C3II and C5 (at 10® yr). These species have been affected by the 
changes in the reaction rate coefficients of O -I- C2, N -I- NO, O + C2H and O -I- C3II 
respectively as they are the most important reaction for each of these species. Among 
these species, the sensitivity of the C5 abundance is an example of the combined effect 
of different rate coefficients since we did not find a unique "key" reaction as done in 
Table|4] The impact of the update in the rate coefficients for dissociative recombination 
is the smallest. 
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Fig. 6 Comparison between model and observations (see text for details). Solid lines are for 
L134N whereas dashed lines are for TMCl. Black lines show the results using the osu-09-2008 
network whereas gray lines are used for the updated network. 
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Table 9 Molecules observed in L134N and TMCl not reproduced by the models at the best 
ages indicated in the table. 



L134N TMCl 



osu-09-2008 


update 


osu-09-2008 


update 


3 X 10^ yr 


2 X 10^ yr 


1 X 10" yr 


5 X 10^ yr 


NO 


CN 


CN 


CN 


HaS 


HaS 


OH 


OH 


CH3OH 


CH3OH 


CH3OH 


C2H 


CH3CHO 


CH3CHO 


CH3CHO 


CH3OH 


HC7N 


HC7N 


C3H3N 


CH3CHO 


HCO+ 


HCO+ 


C3H4 


C3H3N 


HCS+ 


HCS+ 


HCS+ 


C3H4 


SiO 


SiO 


HCNH+ 


HCS+ 






SiO 


SiO 






CH3C3N 


HC9N 






HC9N 





5.3.2 Agreement with observations 

We compared observed and modeled abundances for the dense clouds L134N and 



TMCl, which are normally used as references for these studies (Garrod et al. 20071 



The observed abundances - composed of 42 species in L134N and 53 species in TMCl 



are given in Tables 3 and 4 of Garrod et al. (20071. Upper limits on SiO abu ndance 



are considered in addition (2.6 x 10"^^ for L134N and 2.4 x 10"^^ for TMCl, Ziurys 



et al. 1 1989 1. Since the abundances have been published without estimation of the ob- 
servational error, we assume an error of a factor of three for all molecules. To obtain the 
error bars on the theoretical abundances, we used Monte Carlo simulations to study 



the propagation of uncertainties in reaction rate coefficients (Wakelam et al. 20051. 
The error in the theoretical abundance {2a) is then defined as the envelope containing 
95.4% of the abundance distribution. The agreement between observed and modeled 
abundances is then achieved when an overlap exists between these two quantities tak- 
ing into account both types of error bars. Using this definition of the agreement, we 
computed three parameters as a function of time. The first parameter, labeled F, is 
the percentage of molecules that is in agreement. The second parameter is the dis- 
tance of disagreement D (in log units, similar to Dnorm eq. 26 but not normalized) 
for the species where we do not have agreement. D is the sum, over all such species, 
of the difference between modeled and observed closest error bars. For example, as- 
suming that the observed abundance of a species is defined by Xq^s and Xobs,max 
and the modeled one by X^jj^tj j^j^ and X^-^qj ^^^^ at time t, if the abundance of a 
species is over-estimated by the model, the distance of disagreement for this species 
is |log(Xuiod,min) ~ log(Xobs.max) I ■ Finally, the last parameter allows us to take into 
account both the agreement and the disagreement with the model by multiplying F(t) 
with min(_D)/D(t), with min(_D) the minimum of D over all times. This last parameter 
tends to sharpen the times at which agreement is best. 

We computed the three parameters for four different models. For both clouds, we 
used the same physical parameters: a temperature of 10 K, an H2 density of 10"* cm^'^, 
an Ay of 10 and a ("/^^ of 1-3 x 10"^'' s~^. All species are initially in atomic form, 
except for hydrogen, which is assumed to be already completely molecular. Different 
elemental abundances were chosen for the two clouds. For L134N, we used the typical 
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low-metal elemental abundances (see Table 1 in |Wakelam and Herbst|2008[ ). For TMCl, 
the observation of large quantities of C-rich species seems to indicate that the C/0 
elemental ratio is larger than unity. As in previous studies, we lowered the elemental 
abundance of oxygen to obtain C/0 equal to 1.4. For both clouds, two networks were 
tested: osu-09-2008 and the updated network described in the previous section. Fig. [6] 
shows the three parameters (F, D and F*min(_D)/_D) as a function of time for the two 
clouds and each network. The best agreement for L134N is not changed significantly. 
The best time is still around (2 — 3) x 10^ yr using both networks. For TMCl, we have 
an agreement only slightly better: 79% of the observed species are reproduced with 
the updated network against 77% with osu-09-2008. Specifically, more CH3C3N and 
HCNH"^ are now produced, moving their predicted abundances closer to the observed 
values, whereas C2H is now underproduced by the model. The best times are also 
slightly shifted to earlier times: 4.7 x 10^ — 6.3 x 10^ yr using the updated network 
and 5.7 x 10^ - 1.5 x 10*^ yr with osu-09-2008 although the agreement stays high until 
steady-state. Species not reproduced by the different models are listed in Table [9] for 
both clouds. 



5.3.3 The C + H2 association reaction 

The radiative association reaction C -|- H2 — >■ CH2 -|- hv has been found to be the 
most important one for dense cloud chemistry in the 10^~^ yr era, since it contributes 
significantly to the production of CH2, one of the first carbon-containing molecules. 
The rate coefficient for this react ion recommended in the OSU and UDfA databases is 
lO"'^^ cm'^ and comes from 



Prasad and Huntress 



(1980a I. An estimate based on 



the supposition that radiative stabilisation involves only infrared transitions between 
the vibrational levels of CH2 (see the datasheet in KIDA for more details) yields a value 
of 6 X 10~^^ cm'^ s~^ at 10 K. However, it is possible that the internal levels close to 
the dissociation limit of CH2 are of mixed electronic character allowing an increase 
in the spontaneous emission rate from these levels. Therefore, we have carried out 
calculations with the value of this rate coefficient varied between 10~^^ and 10~^^ cm'^ 
s~^ in order to examine the sensitivity of the results to this reaction. Note that more 
detailed calculations are needed for this reaction in order to have a better estimate. 
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Fig. 8 H2CS and SO abundances as a function of time for cold dense clouds. Different curves 
were obtained with different values of tlie C"'' + S — > S+ + C reaction rate coefficient: 1.5x 10~^ 
(osu-09-2008), f X 10-1°, 5 x 10"" and 1 x 10'" cm^ s'^ 



In Fig. [7j we show the CH2 and HCN abundances as a function of time calculated 
by the updated model when different values of the rate coefficient for formation of CH2 
are employed. HCN is an example of molecules observed in dense clouds and its second 
main production reaction at 10^ yr is N + CH2. The results are shown for five values of 
the C + H2 CH2 + hv rate coefficient: 10"^^, 5 x 10"^^, 10"^^, 5 x 10"^* and lO^^"^ 
cm'^ s~^. The abundances of the species are much more sensitive to the rate coefficient 
when it is larger than 10""'^^. For most molecules, the abundances calculated by the 
model are smaller between 10^ and 10^ yr for smaller rate coef ficients. To reprodu ce 



the observed abundance of HCN in TMC-1, which is 2 x 10"^ (lOhishi et al.||l992|), a 



rate coefficient larger than 10~^^ cm"^ is better. It is however risky to constrain 
rate coefficients by comparing the modeled abundances with the ones observed in the 
interstellar medium since the differences that are found when a single rate coefficient 
is changed will depend on the values of the other rate coefficients and of the reactions 
currently present in the network. The change in other rate coefficients or the addition 
or removal of reactions from the network may change the "best" rate coefficient for C 
+ H2 — >■ CH2 + hv. Detailed quantum calculations would be required to know at least 
if the rate coefficient of this reaction is larger than 10^^^ cm'^s^^ and such detailed 
calculations are beyond the scope of this article. 

5.3.4 The + S charge exchange reaction 

With an accurate quantum chemical treatment for the potential surface, followed by 
semi-classical dynamics, the rate coefficient for the C*" + S charge transfer reaction 
has been calculated to be 5.7 x 10"^^ cm^s~^ at 1000 K and 2 x 10"^^ at 500 K 
( [Bacchus-Montabonel and Talbi|2008[ ). These calculations at high temperature tend to 
show that the 1.5 x 10~^ cm'^s"^ value used in the astrochemical models for the 10-100 
K temperature range is overestimated. A definitive conclusion needs a treatment that 
takes into account the quantum effects at very low temperature, which could increase, 
probably not drastically, the C*" + S rate coefficient. In the meantime, we suggest, for 
temperatures between 10 and 100 K, upper and lower limits for the rate coefficient of 
respectively 10"^" and 10~^^ cm^ s~^. 

In Fig. [Sj we show the abundance of H2CS and SO as functions of time calculated 
for four values of the rate coefficient: 1.5 x 10"^ (osu-09-2008), 1 x 10"^°, 5 x 10"" and 
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Fig. 9 H2CS and SO abundances as a function of time for cold dense clouds. Different curves 
were obtained with different values of two rate coefficients: (1) the rate coefficients of C + H2 
and C+ -I- S are 10"" and 1.5 X 10"^ respectively, (2) C + H2 is lO'^'^ and C+ -I- S is 10" 
(3) C -I- H2 is 10-18 and C+ -I- S is 1.5 X IQ-^ and (4) C -I- H2 is IQ-l^ and C+ -I- S is 10" 
The units of the rate coefficient are cm^ s" 1. 



-10 
-10' 



1 X 10~^^ cm^s~^. The sensitivity does not appear to be very strong. If we, however, 
also vary the rate coefficient of C + H2 , we then have a significant combined effect, 
warned about in the previous section. This is shown in Fig. [9] where the abundances 
of the species have been computed with four sets of rate coefficients (in units of cm'^ 
s"^) for C + H2 and C+ -I- S, respectively: (1) those in osu-09-2008, i.e. 10"^'^ and 
1.5 X 10"^, 2) 10"^'^ and 10^^°, 3) 10"^^ and 1.5 x 10"^ and 4) 10"^^ and 10"^°. 



6 IRC -1-10216 chemistry 



In addition to the chemistry of cold cores, the chemistry of certain circumstellar re- 
gions is also of interest. IRC -1-10216 is a relatively nearby carbon-rich (elemental 
[C]/[0] > 1) asymptotic giant branch (AGB) star which is undergoing a mass loss 
process that produces an extended circumstellar envelope. It is also one of the richest 
molecular sources in the sky with more than 70 molecules discovered to date in its 
circumstellar gas. The hot and dense regions near the stellar surface contain stable 
carbon-bearing molecules such as C2H2, HCN, CS, ... with abundances determined by 
thermodynamic rather than kinetic considerations. As the gas expands it starts to be 
exposed to the interstellar ultraviolet radiation field so that molecules are photodis- 
sociated, and the resulting radicals participate in rapid gas phase chemical reactions 
forming new molecules. The fact that the chemistry of this source is relatively 
simple in terms of geometry (it may be modelled as an uniform spherical 
outflow) makes IRC -|-10216 an ideal target on which to apply and test gas 

19931 



phase chemical models (Cherchneff et al 



Doty and Leung||1998 Millar et al. 2000[ |Agundez and Cernicharo 



Millar and Herbst 



1994 



2006 



Cordiner and Millar "2009). The reader may notice however that uncertain- 



ties in the initial chemical composition (especially the initial abundance of 
C2H2) and UV field may also affect the computed species abundances as 



already noticed by ( jMillar and Herbst 1994 1 
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Fig. 10 Abundance of C4H in IRC +10216 as a function of radius. Left: as calculated with the 
1000 chemical models. Right: mean abundance (solid line) and la standard deviation (dashed 
lines). 



6.1 Chemical model and sensitivity analysis 

Here we report a sensitivity analysis in parallel to the one for cold cores. We computed 
the chemical composition of the gas as it expands from a radius of 10^^ cm to one 
of 10^* cm, while the gas density decreases as the inverse square of radius and the 
temperature as a power law of it with the exponent depending on the distance from 
the central star ( Agundez and Cernicharo|2006 1. 



We utilized the udfaOG chemical network ( Woodall et al.|2007 l with the rate coefE 



cient uncertainties compiled in this database. Because IRC +10216 harbors the largest 
variety of molecular anions yet observed and because of the impact of reactions with 
anions for the chemistry of this source ( Cordiner and Millar||2009 1, the chemical net- 



work was enhanced by a set of reactions involving anions taken from the osu _01_200^ 
database, for which we assumed an uncertainty in the rate coefficients of a factor of 2. 
The photodissociation rate of CO was computed including self-shielding ac- 



cording to the treatment of Mamon et al. ( 1988| ), and the error was assumed 



to be low (a factor 1.25). The sensitivity analysis was carried out by running 1000 
different chemical models in which the reaction rate constants are randomly varied 
within their uncertainties using a lognormal distribution (| Wakelam et al.|2005 |. 



As an example, in Fig.[TO]we show as a function of radius the abundance distribution 
as well as the mean abundance and la deviation of C4H, a species which is formed in 
the outer envelope. The uncertainty in the abundance increases with radius, or with 
time (in an expanding envelope a given radial position may be interpreted as a certain 
temporal instant), which occurs because the dispersion in the abundances due to the 
rate coefficient uncertainties gets larger. 

In the case of a circumstellar envelope such as IRC +10216, it is more appropriate 
to discuss the uncertainty propagation of the reaction rate coefficients in terms of the 
total column densities rather than the abundance at a given radius. In Fig. [lT]we plot 
the total column densities (twice the radial column densities from 10^^ to 10^^ cm) 
calculated with the chemical model for various carbon-chain molecules, with their la 
error due to the uncertainties in the reaction rate coefficients, and compare them with 



^ See http://www.physics.ohio-state.edu/~eric/research.html 
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Fig. 11 Calculated (squares) and observed (circles) column densities of various carbon-chain 
molecules in IRC +10216. 



Table 10 List of key reactions for the chemistry of IRC +10216. 



Reaction i 




FuilC 




Affected species j 


N2 + hi/ N + N 




10 


12.9 


N, N2, HC2N, C3N, C3N-, C2N, SiNH, ... 


C2H2 + hi/ -s> C2H + H 




2 


12.7 


C2H2, C4H2, C3H, C8H2, HC3N, C-C3H2, ... 


C2 + hi/ -s> C^ + e- 




2 


12.4 


C2H+, C2H+, C+, C2H+, H, C2H+, ... 


C2H + hi^ ^ C2 + H 




10 


10.5 


C2H, C3H3, C2, C+, C4H3, C2S, ... 


C + hi/ ^ C+ + c- 




2 


9.3 


C+, e", C, HCP, SiS, H3O+, HCCO+, P, ... 


H2 + cosmic ray — - 


- e 


2 


8.0 


H+, N2H+, HNC, NH+, HCO+, HCNH+, ... 


HCN + hi/ -> CN + H 




1.5 


7.8 


HCN, HC3N, CNC+, C2H3CN, SiH3, HC5N, ... 


C2H + C4H2 CgHa + 


H 


2 


7.1 


C6H2, C7H+, C7H, C4H2, CfiH, ... 


CS + hi^ -s> C + S 




2 


4.3 


CS, S, HC4S+, C4S, H2CS 


Si + hi/ Si+ + e- 




2 


3.7 


Si+, Si, SiNH+, SiN+, SiC3, ... 


C2H + C6H2 ^ CgHa + 


H 


2 


3.4 


C8H2, HC9N, CeHa, C7H, CeH, HC7N, ... 


+ C2 -5- CO + C 




2 


3.0 


C2H+, C2H+, C+, H, C2, C4H+, ... 


C+ + SiS -s> SiS+ + C 




2 


3.0 


SH, SiS+, SO, S2, S0+, SiS, ... 


NH3 + hi/ ^ NH2 + H 




1.5 


2.7 


NH3, NH2, H2CN+, HNO, NH2CN 


C4H2 + hi/ -5- C4H+ + e 




2 


2.6 


C4H+, C6H+, CH3C4H+, C8H+, H3C9N+, C4S 


CN + SiH4 -s> HCN + SiHg 


1.5 


2.5 


SiH4, SiHg, H2SiO, SiCHa, H2SiO+ 


S + hi/ ^ S+ + e- 




2 


2.4 


C4S+, HC2S+, S+, C3S+, NS+ 


CN + C4H2 HC5N + 


H 


2 


2.4 


HC5N, C5N, H2C5N+, CH3C5NH+, CH3C5N 


H2 + SiC2H+ ^ SiC2H+ 


+ hu 


2 


2.3 


SiC2H+, SiC2H2, SiCaH, SiC3H+, SiC3H 



the observed column densities. We note that the observed column densities for long 
carbon-chains are noticeably lower than the calculated ones within their respective 
errors, something which was also found to occur in previous chemical models of IRC 
+10216 ( Millar et al. 2000 1 and which is probably due to incompleteness in the chemical 
network or to wrong reaction rates of key reactions involving these species. The ob- 
served column densities obviously have errors but are less likely to be the cause of large 
disagreements unless some fundamental molecular parameters such as the permanent 
dipole moment are severely wrong. 



6.2 Most critical reactions 



We have used the Pearson correlation coefficients Pj (defined in an analogous manner 
to eq. (25l) with column densities instead of abundances. Table 10 collects the most 
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critical reactions for the chemistry of IRC +10216, as defined by the largest Pj 

where the sum is restricted to species j with a column density in excess of 10^*^ cm~^ 
and with Pj > 0.3. Ifere we have used the second method discussed in section 5.1, i.e. 
the sum of Pj for each reaction. The list of reactions is dominated to a large extent by 
photoprocesses, since the chemistry of the outer envelope of IRC +10216 is very similar 
to that found in photon dominated regions, except that it is a carbon-rich chemistry. 
Also critical are some neutral-neutral reactions as well as the cosmic ray ionization 
of H2, and to a lesser extent ion-neutral reactions. Reactions involving anions are not 
very critical compared with the types of reactions mentioned above, according to the 
adopted criteria. Among them the electron radiative attachment to large radicals to 
form the corresponding anions are the most critical ones, affecting mainly the column 
densities of both the radical and the resulting anion. The CO photo-dissociation 
rate is not very important in our case because of shelf-shielding effects and 
its low uncertainty (F=1.25). 

The types of reactions critical in the circumstellar gas of IRC +10216 are markedly 
different from those for dark clouds such as TMC-1, where ion- neutral reactions dom- 
inate the chemistry and photoprocesses are inhibited by the large visual extinction of 
the cloud. From Table [lO] it seems clear that a critical review of the photoprocesses, 
currently evaluated as a function of the visual extinction through parametric expres- 
sions, is the best initial step in order to reduce the global uncertainty in the chemical 
composition of IRC +10216 as calculated by chemical modelling. 



7 Surface chemistry in models 

Interstellar chemical models that contain both gas-phase and grain-surface chemistry 
are known as gas-grain models. In the standard approach, reactions on grain surfaces 
are treated by rate equations, although the rate coefficients are different from their gas- 
phase counterparts ( Hasegawa et al.|1992 l. More advanced methods utilize stochast 



ic 



approaches, since these are more accurate when the average number of reactive species 
on an individual grain approaches unity, and handle irregular surfaces more appropri- 
ately, as discussed in Section|2] In this discussion, we emphasize the basic rate equation 
method and its uncertainties. In dealing with small particles, it is useful to define con- 
centration in terms of numbers of a given species per grain. Other possibilities are to 
use the bulk concentration of adsorbates and fractional abundances relative to gas- 
phase hydrogen density. The former is obtained by multiplying the number of species 
per grain by the grain concentration. More common units in surface science include 
numbers of monolayers, and areal concentration (cm~^). 

Consider a dust particle with A'^ surface sites, defined as minima in a potential 
surface where adsorption takes place preferentially. Here, we limit consideration to 
the diffusive or Langmuir-Hinselwood mechanism of reaction in which two molecules 
undergo diffusion until finding one another and reacting. Whether the reactants are 
accreted from the gas, or form on the surface from processes such as photodissociation, 
they are assumed to be thermalized before they react. If there are N{A) molecules of 
species A and N{B) molecules of species B on a grain, and if these two species react 
to form molecule C when both species occupy the same potential minimum, the rate 
of formation of C is given by the equation ( [Hasegawa et al.|1992| [Caselli et al.||1998| . 



dN{C)/dt = (fchop.A + fehop,B)iV(^)Ar(B)/A, 



(27) 
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where fc^op.i, the rate of hopping of species i over the potential barrier between two 
adjacent sites, is itself given by the equation 

fchop.i = '^oexp(-Sba/r). (28) 

Here uq is the so-called attempt frequency (s~^) and E\j is the energy barrier in units 



of K. Equation (271 can be simplified to 



dN{C)/dt = (fcdiff,A + kAis,B)N{A)N{B) = kf^BN{A)N{B), (29) 

where fc^jifj j = k]^o-p,i/^ is the hopping rate over a number of sites equivalent to 
the whole grain. If there is a chemical activation energy in addition to the diffusion 
barrier, then a factor k < 1 must be included. The simplest expression for is a simple 
Boltzmann factor containing the chemical activation energy, although more complex 
expressions involving the competition between diffusion and chemical reaction can also 
be utilized ( [Herbst and Millar|2008 L 



Equations for the rate of appearance/disappearance of all surface and gas-phase 
species in a gas-grain network must include terms for adsorption of gas-phase species 
onto dust particles, and desorption of surface species into the gas phase, in addition to 
chemical reactions of formation and destruction. The rate of adsorption, in units of the 
number of species A sticking per single grain per unit time is given by the expression 

dN{A)/dt^Sp^{T)vj^o^,[Al (30) 

where Sa is the sticking coefficient, which is a function of temperature, vp^ is the 
velocity of species A in the gas-phase, crgr is the cross section of the grain, and [A] is 
the concentration of species A in the gas-phase. The equation for desorption depends 
upon the mechanism; for thermal evaporation, which is normally unimportant in the 



cold interstellar medium but increases in importance at higher temperatures ( Collings 
et al.|2004| ) 



it is 

dN{A)/dt = -fecvap,AA^(^), (31) 

where the rate coefficient is the same as that for hopping except that the barrier is 
replaced by the desorption energy Ey^, which is generally considerably larger than E\^: 
fccvap.A = exp (-Ed/T). 

The important parameters in these physical processes are the attempt frequen- 
cies, the energy barriers, and the desorption energies. These parameters are also used 
in stochastic approaches. They are more or less obtainable by a variety of methods 
in surface science, some of which are mentioned below, although estimates are often 
utilized, especially for the diffusion barriers. The result is that calculated abundances 
of surface species are themselves often estimates to an order of magnitude. It would 
appear that the uncertainty and sensitivity analyses discussed for gas-phase reactions 
here can be used for surface reactions as well, but there are many additional com- 
plications including other mechanisms for reaction than the diffusive one, and rough 
surfaces instead of the smooth ones with single values for energy barriers and desorp- 
tion energies. As more is learned about the chemistry on the surfaces of small particles, 
it may well be that uncertainty and sensitivity analyses will be useful. 

Stochastic methods used in place of the rate equation technique can be subdivided 
into macroscopic and microscopic (kinetic Monte Carlo) methods. These are considered 
below in the discussion of H2 formation. 
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8 Surface Chemistry: Experiment 

As an introduction, it must be mentioned that, unlike the situation in the gas phase, 
the results of laboratory experiments on surface reactions are not as easily converted to 
rates under interstellar conditions for a variety of reasons. To study surface chemistry, 
especially at low temperatures, various research groups have performed a variety of 
experiments on cold surfaces such as olivine, amorphous carbon, and assorted ices. In 
general, the experiments are not carried out on small particles. The best characterized 
reaction is the formation of hydrogen molecules from hydrogen atoms ( |Katz et al.| 



19991. The complexities of this process are discussed below. 

In addition to H2 formation, surface reactions in cold cores are thought to produce 
simple ices. Such ices, detected by comparison of astronomical infrared observations 
with laboratory data, are comprised primarily of water, but also contain CO, CO2 and 
several other species, such as CH3OH, NH3, HCOOH, and CH4 in different mixing 



ratios (Ehrenfreund and Schutte 20001. The formation of more complex species on 
surfaces, via a variety of mechanisms, is also likely. In cold sources, reactions between 
atoms that accrete onto the surface and ices may lead eventually to complex species 



(Ehrenfreund and Schutte 2000 Charnley and Rodgors 2005). It is now also thought 



that terrestrially common organic species, such as dimethyl ether and methyl formate, 
found in the gas phase of hot cores, are formed by surface processes involving the 
combination of radicals, whic h are produced pho tochemically from the simpler ices, 
especially on ices of methanol ( Oberg et al. |2009a| . This process is followed by desorp- 



tion as star formation converts a cold core into a hot one ( |Garrod et aL]|2008a| ). The 
experimental study of analogs of these icy grains in space is the topic of a rather new 
research field in laboratory solid-state astrochemistry. 



8.1 H2 Formation on Surfaces 

Molecular hydrogen is the most abundant molecule in the ISM. It is a precursor for 
the formation of more complex molecules through protonated molecular hydrogen, , 
serves as an important coolant, and shields itself and other molecules from the harsh 
interstellar radiation field. This shielding arises from the unusual two-step photodis- 
sociation mechanism of H2 in the Lyman and Werner transitions from 912-1100 A, 
which is discrete rather than continuous, so that specific small wavelength regions are 
depleted strongly. Under most circumstances, the most efficient formation route of H2 
in the gas phase is through three-body reactions where the third body can take some 
of the 4.5 eV of excess energy that is released in the reaction. In the dilute conditions 
of the ISM, three body collisions do not occur and these reactions are not possible. 
Already in 1949 it had been suggested that surface processes on dust particles could 

. The grain surface would in this case 
act as the third body. 

The surface formation of molecular hydrogen can roughly be divided into two 
regimes. At low temperatures, the hydrogen atoms are weakly bound, or physisorbed, 
to the grain surfaces and can still easily move on the surface. Once two atoms meet, 
they will form molecular hydrogen through a barrierless reaction. If the surface tem- 
perature becomes too high, the residence time of atoms on the surface becomes so short 
that, at the low fluxes in the ISM, the chance of two atoms meeting becomes negli- 
gible. At these higher temperatures, H2 formation has therefore to proceed through 



play a role in H2 formation (van de Hulst|p49l 
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chemisorbed hydrogen atoms. Observations show indeed that molecular hydrogen is 
also formed in warmer areas (> 20 K) such as PDRs and post-shock regions. Here we 
will make a clear distinction between the two mechanisms and we will discuss the main 
uncertainties related to both cases. 



8.1.1 Low temperature 



The H2 formation process from physisorbed atoms has been intensely studied over the 
past years. These studies have focused on three different surfaces: silicates, carbona- 
ceous material and water ice. Experimental studies are geared towards obtaining bind- 



ing energies ( 


Pirronello et al.||1997a|b 


1999 


Hornekaer et al.||2003 2005 Perets et al. 


2005 Dulieu et al.||2005 


Matar et al. 


2008 


Pillion et al.||2009 


1 or excitation energies 


of the desorbing H2 molecules ( 


Perry and Price 


2003 


Congiu et al. 20091. Theoret- 



ical studies have used the binding energies as input parameter and have focused on 



the influence of surface roughness and amorphicity ( Chang et al.|2005 Cuppen et al. 
2006 Cuppen and Herbst|2005" I and the so-called accretion limit (Charnley 1998 2001 
Stantcheva ct al. 2002| [Vasyunin et al.|2009T |Biham et al.||2001| [Green et al.||2001| ) 011 

the formation efficiency. 

Experimentally, molecular hydrogen formation is studied using mass spectrometry 
coupled with temperature-programmed desorption (TPD). In an ultra-high- vacuum 
set-up, hydrogen atoms are deposited on a cold interstellar grain analogue of one of 
the above-mentioned materials and the desorption of new molecules as a function of 
temperature is measured. Usually beams of partially dissociated H[2 and D2 are used 
and HD is observed, thereby avoiding the problems associated with the presence of 
undissociated molecules in the incident beams. These beams are either produced in a 
micro-wave or a thermal cracking source. The atoms are generally coUisionally cooled 
before they reach the surface where they can react. Binding energies can be obtained 
by fitting a rate equation model to the experimental data where the desorption and 
diffusion energies are used as fitting parameters. The procedure results in different 



types of uncertainties. In their work, Katz et al. ( 1999 1 stated that they were only able 



to extract lower limits for the desorption energy of H atoms from the experimental data 
of [Pirronello et al.| ( [1997a|b| [i999J . Higher values for this energy were found to give 
degenerate results. In general, one could state that the desorption energies of stable 
species like H2 and D2 are relatively straightforward to determine whereas desorption 
energies of unstable species like H and D are much harder to determine. Unfortunately, 
the latter are especially important in the determination of the temperature window in 
which efficient formation of H2 can occur. 

Another problem with the fitting procedure is that the formation mechanism is 
implicitly assumed in the model used for fitting. To use the results from Katz et al.| 



( 1999 1 again as an example, Cazaux and Tielens ( 2004 1 and Cuppen and Herbst ( 2005 1 



were able to fit the same experimental results using different models. Cazaux and| 
[Tielens[ ( [2004[ ) applied a binary model with both physisorption and chemisorption 
binding sites, as well as quantum mechanical tunneling and thermal hopping between 
sites. Cuppen and Herbst ( 2005 1 introduced surface roughness in their model. Another 



crucial ingredient to H2 formation models is the diffusion barrier, which was recently 
determined to be 22 ± 2 meV for D atoms on porous amorphous water ( jMatar et al.[ 
2008 1 in agreement with earlier experimental ( Hornekar et al.][2003 l and theoretical 
results ( Buch and Czerminski[[l991 1. The determination of diffusion barriers remains 
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experimentally challenging, however, and is one of the largest causes of uncertainty in 
the modeling of surface chemistry in general. 

The experimentally determined desorption rates show a rather large spread. This 
spread is not only due to experimental uncertainties but is an intrinsic property of the 
interstellar grains. [Hornekar et al.| ( |2003[ |2005[ ) and |Fillion et al.| ( [2009[ ) have shown 
the importance of surface roughness and porosity for H2 desorption from amorphous 
solid water ice. More porosity or roughness leads to a wider distribution of desorption 
energies with a higher mean value. The exact preparation of the water ice substrate 
is key to the amount of porosity. [Perets and Biham| ( |2006[ ) showed a similar effect 
occurs for silicates. Computer simulations of molecular hydrogen formation under in- 
terstellar conditions, i.e., using much lower fluxes than under laboratory fluxes, had 
already shown that an increase in surface roughness allows H2 formation at higher tem- 
peratures ( [Chang et al.|[2005l [Cuppen et aL||2006| [Cuppen and Herbst||2005| . Surface 
irregularities in the form of protrusions, pores or amorphicity can act as sinks which in- 
crease the residence time of atoms on the surface whereas atoms on the smooth terraces 
can easily move around and reach these sinks. A second large uncertainty in surface 
chemistry is therefore the uncertainty in the exact surface grain structure. It appears 
that the difference in binding energies between the different materials is smaller than 
the difference due to the differences in surface structure. 

Unlike the surfaces used in the laboratory, interstellar dust grains are small and 
generally have a small number of H atoms attached to them. The average number of H 



atoms per grain can under certain interstellar conditions be unity or less ( Tielens and 



|Hagen|1982[ [Biham et al.|200T| [Lipshtat et al.|2004a[ ) . Since clearly a minimum of two 
atoms needs to be present for H2 to form, deterministic models like the standard rate 
equations can break down and one has to resort to stochastic methods. Generally two 
types of models are applied to replace the rate equation methods: macroscopic Monte 
Carlo methods ( |Charnley|1998l[200T||Stantcheva et al.|2002 !'Vasyun in et al.|2009|) and 



master equation treatments (Biham et al 



2001 



Green et al. 2001, Stantcheva et al 



2002 1 . The latter can be easily linked to rate equations which handle the gas phase 



chemistry in the most straightforward way. The main disadvantage of this technique is 
that, as the number of species expands, the equations blow up rapidly. Algorithms have 



been suggested to make the effect less dramatic (Lipshtat and Biham 2004b Barzel 



and Biham,2007| ). Apart from the stochastic methods, some modifications to the rate 
equations in a semi-empirical manner have been put forward to mimic some of the 



results of these methods ( Caselli et al. 1998 2002 Garrod 2008b I. The above-mentioned 



roughness can be included via microscopic or kinetic Monte Carlo methods ( Herbst and 



Cuppen|2006| ), and appears to extend the regime in which deterministic methods apply, 
since it increases the surface abundances. Whether simple rate equations will be able 
to capture this roughness is not yet clear (Cuppen et al. 2010a, submitted). 



8.1.2 High temperature 



Molecular hydrogen formation at higher temperatures is almost uniquely studied on 
graphite. This system acts as a model system for both carbonaceous grains and poly- 
cyclic aromatic hydrocarbons (PAHs). Hydrogen atoms can form a chemical bond with 
graphite after they overcome a chemical activation barrier of 0.15-0.2 eV ( |Hornekaer| 
et al.|2006b Jeloaica|1999 1. Chemisorbed atoms therefore mostly occur at higher tem- 



perature. 
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Experimental (|Hornekffir et al.||2006a 


Andree et al. 


2006 


1 and theoretical (Ferro 


2003 Miura et al.||2003 


Hornekffir et al. 


2006a|b 


Rougeau et al. 2006 1 studies show 



that hydrogen atoms preferentially adsorb in clusters on graphite. Density Functional 
Theoy (DFT) calculations show that not only the sticking into dimer configurations is 
more favorable, but also that the binding energy of these dimers is higher than for two 
individual monomers. More complex structures consisting of three or more hydrogen 
atoms are found to be favorable as well (Hornekaer et al. '2006b Casolo et al.|2009 1 . The 



complex desorption behavior of H2 from graphite as observed by Zecho et al. 
was attributed to the presence of dimers and at high coverage to tetramers 
|et al.|2009| ) 



(2002a I 



Gavardi 



Hydrogen abstraction from graphite has been studied experimentally ( Zecho|2002b 



Ze^aetalj2002a \ and by DFT and Quantum Wave Packet calculations ( Farebrother 
et al.|2000||Meijer et al.|2001||Slia[2002l|Sha et al.|2002[[M"orisset et al.|2004[|Martinazzo 



and Tantardini 2006 \ . Since the barrier for diffusion is higher than the desorption 
barrier for chemisorbed atoms, the formation mechanism is often referred to as an 
Eley-Rideal reaction where an incoming hydrogen atom reacts with another H atom 
that is chemically bound to the surface. Computational limitations do not, however, 
allow a simulation of the full complexity of the H-graphite system. Simplifications 
could be problematic. Generally either none or only a single C atom on the simulated 
graphite surface is allowed to relax during the reaction, which could artificially close low 
barrier reaction channels, such as sticking of hydrogen atoms into specific configurations 
( Hornekaer et al.]|2006b Rougeau et al. 2006 ). Cuppon and Hornekaer (2008j) showed, 
in a kinetic Monte Carlo study which did allow for different mechanisms, that H2 is 
more likely formed through a combination of fast diffusing physisorbed atoms ( |Bonfanti| 
et al. 20071 and chemisorbed atoms. The fast moving physisorbed atoms can find a 
chemisorbed atom to form a dimer with. Part of the excess energy released in the 
binding process can be used to overcome the barrier for H2 formation and desorption 
from the formed dimer. 

Under interstellar conditions, we expect some H2 to form upon H atom sticking 
through this dimer-mediated mechanism. If the cloud is too cold for the remaining 
complex structures to desorb, stochastic heating of the carbonaceous grains by UV 
photons can facilitate this process. Detailed simulations under these conditions are in 
progress and indeed the dimer-mediated mechanism proves to be dominant for high 
(gas) temperature (Cuppen et al. 2010b, submitted). 



8.2 Ice chemistry 

The formation of H2 occurs in the submonolayer regime. Molecular hydrogen does not 
build up in layers but easily desorbs. When less volatile species land or form on the 
surface, layers of ice form, which can be rather inhomogeneous in composition. This 
layering gives an extra dimension to the surface chemistry which makes the laboratory 
results harder to interpret and surface modeling less intuitive. 

The new research field of laboratory solid state astrochemistry involves surface 
processes such as desorption, reactions and energetic processing. Again models are 
used to interpret these data. Information on one-step processes like thermal desorption 
and photodesorption is relatively straightforward to extract. However often a sequence 
of processes as well as competing processes make it hard to disentangle the individual 
steps that are used as input for grain surface models. Here we will discuss thermal 
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desorption, surface reactions by means of atom bombardment, and UV photoprocessing 
of ices. 



8.2.1 Thermal desorption 



As in the case of hydrogen, the desorption energies determine the temperature regime 
in which species will be present on the surface of a grain and available for reaction. The 
desorption energies of a wide collectio n of stable species have been determined u sing 
the TPD technique. E xampl es are N2 (|6berg et al.|2005|, CO |6berg et al.|2005|, O2 
(Acharyyaet al.|2007[), H2O ( [Bolina et al.|2005[|Fraser et al.|2001| ), and CH3OH ( [Green 



et al.|2009 K The desorption energies have been mostly determined for the desorption of 



pure ices from different substrates. The difference between the different substrates are 
rather small and become negligible in the multilayer regime. In the multilayer regime, 
the molecules desorb with a (near) zeroth order rate whereas they desorb with a (near) 
first order rate in the monolayer regime. 

Since interstellar ices are not homogeneous the desorption of mixed layers is more 
relevant for astrochemical modeling. However, the introduction of more species in the 
ice makes the desorption process immediately much more complex. First, the desorption 
energy can change depending on its surrounding material. Second, the dominant mantle 
species can prevent other species from desorbing. Collings et al. ( 2004 1 showed, for 
instance, that molecules like CO and CO2 can become trapped in an ice mantle which 
consists predominantly of water ice since the desorption of water occurs at much higher 
temperatures than for CO and CO2. However, at the long timescales available in the 
ISM, some of these trapped species might be able to escape be cause of segregation , 
in which the two main fractions of the mantle slowly separate (Oberg et al. 20091 1. 
This process depends on a large number of parameters including surface temperature, 
ice composition and mixing ratio. Fortunately, it appears that s ome simple ch anges in 

FayoUe 



2009d 



gas-grain codes could treat the main features of this behavior ( Oberg 
et al. in preparation). These changes would include splitting the grain mantle into 
two phases: a bulk and surface phase ( Hasegawa and Herbst|1993 1. Desorption is only 
allowed from the latter phase. 



8.2.2 Surface reactions by atomic bombardment 



The experimental studies of ice reactions can be roughly divided into two groups de- 
pending on the analysis technique and the thickness regime (submonolayer versus mul- 
tilayer). As in H2 formation, the reactants and products can be probed mass spec- 
trometically. The surface is initially exposed to a small quantity of the reactants, after 
which the surface is heated until the products and reactants desorb and are detected 
(thermally programmed desorption or TPD). Different initial exposures and tempera- 
tures can be probed to obtain information on reaction order, etc. The main advantage 
of this technique is the sensitivity, which allows for submonolayer exposures and which 
is able to detect all species by their masses. The method has however four major disad- 
vantages: the products cannot be probed in-situ, i.e., during the atom bombardment, 
additional reactions during the heat-up to desorption cannot be excluded, quantifying 
the desorbing species is not straightforward, and some of the interesting species have 
equal masses. 

A second method is to initially grow an ice of several monolayers and expose this ice 
to the atom beam while Recording Reflection Absorption Infrared Spectra (RAIRS). 
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In this way the reactants and products are probed in-situ at the time and temperature 
that one is interested in, which is the main advantage of this technique. Quantifying 
the formed product is relatively simple, provided that the RAIRS is calibrated with an 
independent method. The main disadvantages are that not all species can be detected 
unambiguously, e.g. because of spectral overlap, and that the sensitivity is less than 
with the previous technique. So far this technique has only been applied to smaller 
species, with ethanol the largest product formed (I Bisschop et al.|[2007 K In the recent 
years, most laboratory reaction studies have been focused on the formation of two 
molecules: methanol and water. We discuss both briefly below. 

Let us start with the formation of surface methanol, which has long been thought 
by astrochemists to proceed by the hydrogenation of CO ice via successive reactions 
with atomic hydrogen ( Charnley et al.|1997 |: 

CO HCO H2CO ^ H2COH ^ CH3OH. (32) 

Laboratory data have become available confirming this sequence of reactions in several 



laboratories using RAIRS as the main analytical technique (Watanabe and Kouchi 



2002 Fuchs et al. 20091. Generally, pure CO ice is grown under ultra high vacuum 
conditions (< 10"^" mbar). The ice is subsequently bombarded by H-atoms, which 
are produced by dissociating molecular hydrogen, II2, using a microwave or thermal 
cracker. The formation of H2CO and CH3OH was detected by both RAIRS and TPD. 
The first and third reaction step possess barriers to reaction. Since the rate of formation 
of molecules depends on a sequence of events: deposition of H atoms, diffusion and then 
reaction, these barriers can only be determined using a detailed model. Moreover, the 
formation of H2CO and CII3OH is in competition with H atom desorption and H2 
formation. A product yield of a few monolayers was obtained, which increases with 
temperature, below the desorption temperature. This indicates that the H atoms are 
able to penetrate a few monolayers of the CO ice and that this penetration depth 
increases with temperature. A microscopic kinetic Monte Carlo algorithm was fitted 
to a large dataset of H2CO and CH3OH production and CO loss data as functions of 



temperature and H-atom fluence (Fuchs et al. 2009 1. The reaction barriers were used as 



fitting parameters. The uncertainties in the obtained chemical barriers are mostly due 
to uncertainties in diffusion barriers and sticking fractions and experimental conditions 
such as temperature and fluxes. Using these barriers, the methanol production can be 
simulated for interstellar conditions ( Cuppen et al.|[2009 |. This study clearly showed 
that under interstellar conditions different processes become dominant as compared to 
the laboratory. Under laboratory conditions H2 formation is dominant whereas under 
interstellar conditions CO hydrogenation is dominant. The study further showed that 
it is crucial to include layering and consider accessibility in surface models in order 
to correctly simulate surface chemistry. Again a model that distinguishes between the 



surface layer and internal layers could help to include this behavior (Hasegawa and 
Herbst|1993'l ) 



Water ice can be formed on grain surfaces through three reaction channels: hy- 
drogenation of atomic oxygen, molecular oxygen and ozone ( [Tielens and Hagen '1982]). 
The laboratory of the Universite Cergy-Pontoise has studied all three channels using 
a Quadrupole Mass Spectrometer as their main analysis technique ( Dulieu et al.|2007 



Matar et al. 2008 Mokrane et al. 20091. They confirmed that the hydrogenation of 



these three species indeed results in the formation of water, but quantitative data are 



not available. The hydrogenation of solid O2 is the best studied channel (Miyauchi 



et al. 2008 loppolo et al. 2008 Oba et al. 20091. loppolo et al. (20081 investigated 
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this reaction for a large range of temperatures. They found very surprising results. Un- 
like the CO hydrogenation study, where the product yield was only a few monolayers; 
O2 hydrogenation appeared to result in a much larger yield which showed a strong 
temperature dependence. This larger yield occurs because of the ability of H atoms to 
diffuse easily through an O2 ice (loppolo et al. in preparation). These results compli- 
cate the layering problem as discussed above, since the shielding of molecules is found 
to be strongly dependent on the species. Furthermore, the experiments showed that 
the reaction network of three channels that had been initially proposed appears to be 
more complex and that several of the channels are linked through additional reactions 
(Cuppen et al. 2010c, submitted). 



.3 UV-photoprocessing of ices 



Other experiments are directed at the chemical and physical processes that occur dur- 
ing and after UV processing of ice samples. The flux of UV photons in the laboratory is 
much greater than the flux in space, so that great ca re must be taken in the i nterpreta- 
tion of these experiments. In the Leiden experiment ( Oberg et al. |2005| |2007[ ) , the ice is 
irradiated by the output of a broadband hydrogen-microwave discharge lamp that sim- 
ulates the interstellar radiation field. It peaks around 121 nm (Ly a) and covers 115-170 
nm (~7.5-lQ.5 eV) with a variable photon flux covering roug hly (1 - 10) X lO" pho- 
tons cm^^s~^ at the substrate surface. The resulting ice behavior is monitored using 
RAIRS to quantify changes in ice constituents and TPD to investigate ice desorption 
and identify the products of ice chemistry. 

One of the major outcomes of these stuides is that photodesorption, a non-thermal 
desorption mechanism, is quite effi cient and thus important in the cold interstellar 

for CO, N2, CO2 and H2O ice photodesorption 



medium (see Oberg et al 



2007 



2009b|c 



efficiencies, respectively). The investigated molecules photodesorb through different 
mechanisms. N2 ice only co-desorbs with other molecules, which is expected since it 
has no dipole-allowed transition within the wavelength range of the lamp, while CO ice 
photodesorbs non-dissociatively. H2O and CO2 are dissociated by the lamp. Following 
dissociation, theoretical simulations predict that a fraction of the fragments will desorb 
directly or recombine in the ice and then desorb (Andersson and van Dishoeck||2008[ ). 
Both desorption pathways are observed experimentally. 

The photodesorption yields are ~10~^ molecules per incident photon for all molecules 
except for N2, which has a yield that is an order of magnitude lower. CO ice photodes- 
orption is a surface process because of its peculiar mechanism. CO2 and H2O are 
observed to photodesorb from deeper within the ice, but desorption is still limited to 
the top 10 layers. It is also found that for methanol the photod esorption yield amou nts 
to around 2x10"'^ desorbed molecules per incident UV photon (Oberg et al. |2009e| . If 
one looks at the photodesorption in more detail, the yield Y does depend to some extent 
on T, the ice temperature, x, the ice thickness an d IjT), the temperatur e dependent 



mean-free-path of the excited molecule in the ice. |6berg et alT] ( |2009b|c| empirically 



determined photodesorption yields Y for CO, CO2 and H2O to be 



ico = 10 ^ (2.7- (T- 15) X 0.17), (33) 

-1/2. 9x 



ic02 (T <AOK)^ 10-^ (1.2(1 - e"^/2-^) + 1.1(1 - e"^/''-^)) , (34) 

" (^(1.3 + 0.032 X r)(l - e"^/'(^))^ . (35) 
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These yields are for multi-layered ices and have absolute uncertainties of ~40-60%; 
the current experimental set-up does not allow for (sub-)monolayer desorption deter- 
minations. Despite their uncertainties, these relatively high photodesorption rates may 
explain, for example, why CO ca n be found in the gas ph ase at temperatures well 



2009b I , in protoplanerary disks 



below its desorption temperature (Oberg et al. ] |2007| 
for instance ( [Hersant et al.|2009| ). 

As photodesorption proceeds, photodissociation also occurs and new spectral fea- 
tures appear that correspond with reaction products ( Oberg et al. |2009a| . Photochem- 
istry in icy grain mantles was suggested as a path to chemical complexity more than 



three decades ago (Hagen et al. 19791 and studied since then (D'Hendecourt et al. 



1986 Dartois et al. 



2005 



Nuevo et al. 20071. New experiments in Leiden show that 



many of the complex organics observed in hot cores indeed can be formed through UV 
irradiation of solid CH3OH and ice mixtures containing methanol ( Oberg et al. |2009e| . 
The UV radiation dissociates methanol into CH2OH, CH3O, CH3, OH, and H. De- 
pending on ice temperature, some or all radicals react to form more complex molecules 
such as C2H5OH and (CH20H)2. During star formation, ices are always exposed to 
low levels of UV irradiation because of cosmic ray interactions with H2. After the pro- 
tostar turns on, the temperature of the cold material rises and the radicals produced 
by the photoprocessing diffuse rapidly enough to react to form more complex species, 
which are thermally released into the gas phase ( [Garrod and Herbst|2006| [Garrod et al. 



2008a I. The laboratory results are currently being fit to a theoretical model that uses 
rate equations; once the basic rate parameters are understood, they can be used to 
refine the model of photo-production of complex organics in hot cores (Oberg, private 
communication) . 



9 Summary and conclusions 

Some of the major projects undertaken in connection with this article are: 

(i) the identification of the 'key' reactions in gas-phase models of the cold cores in 
dense interstellar clouds and of circumstellar envelopes; that is, those reactions that 
influence most strongly the predicted abundances of the major species observed in 
these molecule-rich astronomical objects, and their precision; 

(ii) the careful evaluation of the available information on both the rate coefficients 
for these processes and the uncertainties in these coefficients; 

(iii) an examination of the effects of including the re-evaluated rate coefficients on 
the predicted abundances of the gas-phase model of cold cores. 

In parallel with the preparation of this article, a new database for astrochemistry 
(| http :/ /kida.obs.u-bordeauxl.fr/ ) has been prepared. Reasons for choosing the new 
values of rate coefficients can be found in that database. In running the models with 
the newly evaluated input data, uncertainties in the predicted abundances have been 
estimated by running a large number of model calculations with the rate coefficients 
for each of the reactions in the model chosen randomly within the evaluated uncertain- 
ties. The impact of changes in the rate coefficients has been explored for two types of 
astronomical sources, in which the physical conditions are quite different: cold dense 
clouds and evolved stars. In the latter case, photodissociation processes, which have 
not been evaluated here, are shown to be important. The comparison between pre- 
dicted and observed abundances is slightly improved for the TMC-1 dense cloud by 
including the new rate coefficients. Whatever the current degree of improvement, it is 
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important to include the 'best' values of the rate coefficients in this, and future, model 
calculations. Outstanding amongst the key reactions is the radiative association of C 
atoms with H2 to form CH2. This reaction presents an especially difficult problem to 
both experimentalists and theoreticians, but one that we hope will be tackled in the 
near future. 

In the latter part of our article, we discuss how chemical changes on the surface of 
interstellar grains can be incorporated into so-called gas-grain models. The advantages 
and disadvantages of both rate equation and stochastic methods are considered. In ad- 
dition, the results of recent experiments on surfaces designed to mimic the chemistry 
on both 'bare' surfaces and those 'coated' with ices of different composition are re- 
viewed, with some emphasis placed on the mechanisms for the formation of molecular 
hydrogen. 
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A Appendix A 



Author: Ian Smith (University of Cambridge, UK) 



H3* + 0fP) ^ 0H" + H2 (1) 

^ HzO^ + H (2) 

Thermodynamic Data 

AH°29s(l) = - 57 kJ mol-' Afl°298(2) = - 158 kJ mol"' 

Thermochemical data are taken from ref. 3 

Rate Coefficient Data {k=ki + Icz) 



k/cw? molecule"' s"' 



Reference 



Comments 



Rate Coefficient Measurements 
(8±4)X lO '" 
(12±5)x 10'" 
Theoretical Evaluations 

11.9 X 10-'" (r/Kysoo)-"'"" 

9.4 X 10-'" (77K)/300)-°^'° 



300 

295 ±5 

0-70 
70-500 



1.14 X 10-' (77K)/300)-" '5'^exp(-1.41/7) 5 - 400 



Fehsenfeld, 1976' (a) 

Milligan and McEwan, 2000^ (b) 

Bettens, Hansen and Collins, 1 999" (c) 

Bettens, Hansen and Collins, 1999" (c) 

Klippenstein et al. , 20 1 0' (d) 



Reviews and Evaluations 
(8±4)x 10'" 
(8±4)x 10-'" 



all temperatures UMIST database 
all temperatures OSU website 



Comments 

(a) A very short paper with few experimental 
details. Flowing afterglow technique: 
production of Yl^ described in ref 3, no 
details of how of 0-atom concentration were 
estimated. No measurements of branching 
ratio. 

(b) Experiment employs a SIFT technique to 
generate and select Hs^ ions. O atoms from 
microwave discharge of Oi/He. Reaction of 
CH3^ + O (which it is claimed has a well- 
established rate coefficient) was used to 
calibrate the 0-atom concentrations. Error 



takes into account uncertainties in estimating 
0-atom concentrations. Branching ratios are 
found to be k\l{ki + fe) = 0.7 and ^2/(^:1 + k^) = 
0.3. 

(c) Results of a quasiclassical trajectory study 
on an ab initio potential energy surface/ 
provides a temperature-dependence for k. The 
branching ratios were calculated to be k\l{k\ + 
ki) = 0.94 and ki/ih + ki) = 0.06, independent 
of collision energy. 

(d) Kinetics predicted from a combination of 
transition state theory, trajectory calculations 
and master equation analysis employing 



potential energy surface from high level 
quantum chemical calculations. 

Preferred Values 

;t(298 K) = (1.2 ± 0.3) X 10"' cm^ molecule"' s" 

1 

kilO K) = (1.7 ± 0.4) X 10-" cm' molecule' s"' 
k(T) 1.14 X 10-" (r/K)/300)-° '''^exp(-1.41/7) 
cm^ molecule"' s"' 

ki/(ki + ki) = (0.7 ± 0.2) and + fe) = (0.3 
± 0.2), independent of temperature. 



Reliability 
AlogA: = ±0.15 

Comments on Preferred Values 
Both the OSU and UMIST databases accept 
Fehsenfeld's room temperature value from 
1976,' and (apparently) assume (a) that 
reaction (1) to OH^ + H2 is dominant, and (b) 
that the rate coefficient is independent of 
temperature. Unfortunately, Fehsenfeld gives 
no details of how the 0-atoms were generated 
and their concentration estimated. A follow- 
up paper that was supposed to give these 
details could not be found. 

Milligan and McEwan^ do give details. In 
effect, they measure the rate coefficient 
relative to that for CHjO^ + O. Their 
assessment of errors seems well founded, 
their value of ^(298 K) and the branching 
ratio at 298 K (and the uncertainties in both 
quantities) are accepted. 

Both Fehsenfeld' and Milligan and McEwan^ 
only give room temperature values of the rate 
coefficient. The latter authors point out A(298 
K) is close to the value expected from the 
Langevin model (17 x 10"'° cm' molecule"' 
s"'), on which basis one might expect no 
temperature-dependence. 



Bettens et al.'^ provide a thorough theoretical 
analysis of this reaction. Their findings include: 
(a) that the capture cross-section is greater than 
that given by the Langevin model due to the 
long-range ion-quadrupole forces between 
and 0('P); and (b) that only ca. 80% of the 
'captured' trajectories lead to reaction. They 
assume that reaction occurs only on the lowest 
triplet surface correlating with reactants and 
reduce their calculated value by a factor taking 
account of the spin-orbit states in 0('P,/). 
Together with (a) and (b), the value of A:(298 K) 
obtained is close to that obtained by Milligan 
and McEwan^ and to that estimated by a 
Langevin analysis. 

The calculations of Klippenstein et al.^ are in 
fair agreement with those of Bettens et al.'', 
both in respect of the absolute values of the rate 
coefficients and their T-dependence, and again 
with the limited experimental data. They do not 
provide information on the branching ratio. We 
adopt their recommended T-dependence (in 
recognition of their more sophisticated 
treatment of spin-orbit effects) and the 
branching ratio from ref. 2. 

We adopt the values of A:(298 K) and the 
expression for k(T) from Klippenstein et al. In 
view of the good agreement between the two 
sets of calculations and their agreement with 
the experimental values of k(298 K) the 
uncertainties in the recommendations are quite 
small. 
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